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ABSTRACT 


A  balloon-borne  aerosol  counter  which  samples  and  sizes  individual 
particles  to  an  accuracy  of  +  10  percent  over  an  order  of  magnitude 
size  range  extending  down  to  the  0.1  -  0.2  micron  diameter  regime,  has 
been  successfully  flown.  The  device  is  completely  self-contained  and 
recoverable.  The  data  are  recorded  on  9-channel  magnetic  tape  in  digi¬ 
tal  format  compatible  with  an  IBM  360  computer,  and  are  processed  auto¬ 
matically. 

Stratospheric  aerosol  size  distributions  have  been  obtained  in 
increments  of  one  and  two  kilometers,  up  to  and  including  balloon  float 
altitudes  (above  23  km).  These  represent  the  first  such  direct  infor¬ 
mation  available,  and  are  presented  here  in  preliminary  form. 

A  description  of  the  device,  together  with  a  discussion  of  cali¬ 
bration  and  automatic  data  reduction  procedures,  as  well  as  recommenda¬ 
tions  for  improvements  in  the  latter,  is  given. 
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A  BALLOON-BORNE  AEROSOL  COUNTER 


I.  INTRODUCTION 

This  Final  Report  under  Contract  FI9628-70-C-0265  describes  the 
preparation  for  and  successful  test  flight  of  a  balloon-borne  aerosol 
particle  counter.  The  device  samples  individual  particles  and  measures 
their  sizes  to  an  accuracy  of  +10  percent  over  an  order  of  magnitude 
extending  down  to  the  0.1  to  0.2  micron  diameter  regime.  The  size 
distribution  of  aerosols  at  various  stratospheric  altitudes  during  the 
ascent  and  descent  portions  of  the  balloon  flight,  as  well  as  at  float 
altitude,  was  measured  successfully  on  the  first  balloon  flight  on  28 
October  1970  at  White  Sands  Missile  Range.  Preliminary  results  of  this 
flight,  which  to  our  knowledge  represent  the  first  such  data  obtained  to 
date,  are  presented  here.  This  device  was  designed  and  fabricated  under 
a  previous  program,  Contract  No.  F19628-68-C-0086,  but  had  sustained 
severe  damage  prior  to  its  first  flight  test  as  a  result  of  its  accidental 
exposure  to  temperatures  in  excess  of  200°C  in  an  environmental  test 
chamber.  This  necessitated  the  postponement  of  the  originally  scheduled 
field  exercise,  and  the  initiation  of  an  extensive  program  of  renovation 
and  refurbishment. 

A  description  of  the  aerosol  counter,  as  well  as  a  damage  assest  u-nit 
report  and  recommendations  for  renovation  has  been  presented  in  Reference  1. 
These  recommendations,  which  included  certain  specific  modifications  to 
the  ' '.act  ronics  circuitry  and  especially  to  the  ground  support  system, 

.  V  n  implemented  under  the  present  contract,  and  are  described  sub¬ 
sequently  in  the  present  report. 

Most  of  the  mechanical  and  optical  components  were  repaired  under 
the  prior  contract  in  order  to  facilitate  the  damage  assessment  to  the 
electronics  circuitry.  This  was  necessitated  by  the  unique  input  require¬ 
ments  of  the  several  digital  modules  involved.  Hence,  the  renovation 
work  performed  under  the  present  contract  has  been  devoted  largely  to  the 
electronics.  During  the  course  of  this  activity  several  additional  modi¬ 
fications  and  improvements,  which  had  not  been  contemplated  originally, 
were  incorporated  for  reasons  which  are  discussed  in  greater  detail 
herein. 

For  the  benefit  of  the  reader  who  may  not  have  ready  access  to 
Reference  1,  a  brief  synopsis  of  the  counter  system  is  presented  in 
Section  2  of  this  report. 

A  thorough  discussion  of  all  renovation  activity,  including  calibra¬ 
tion  and  environmental  testing,  is  presented  in  Section  3.  With  regard 
to  calibration,  it  had  been  recognized  that  the  scheme  employed  earlier, 
namely  the  use  of  Dow  Corning  "standard"  polystyrene-latex  spheres, 
though  adequate  in  terms  of  providing  a  suitable  range  of  particle  sizes, 
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was  nevertheless  deficient  in  that  only  one  index  of  refraction  was  available, 
namely  1.58.  Under  the  present  program,  a  technique  was  employed  to  generate 
particles  of  a  known  size  with  index  of  refraction  equal  to  1.50.  This 
rendered  it  possible  to  test  the  capability  of  the  aerosol  counter  to  measure 
particle  sizes  accurately,  (to  within  the  prescribed  10  percent),  over  the 
anticipated  range  of  particle  index  of  refraction  (1.33  to  1.50).  The 
results  obtained  serve  to  strengthen  one's  confidence  in  this  particular 
feature . 

A  description  of  the  flight  test  activity  is  presented  in  Section  4. 

In  addition,  detailed  discussions  of  both  the  rationale  employed  for  the 
automatic  processing  of  scientific  data  as  well  as  a  complete  description 
of  the  manner  in  which  this  and  other  vital  information  is  recorded  on 
tape  are  also  included. 

The  aerosol  counter  was  flown  on  two  successive  balloon  flights  during 
a  two  week  period  in  October  1970  at  White  Sands  Missile  Range.  Excellent 
data  were  obtained  on  the  first  of  these  flights,  while  all  data  were  lost 
on  the  second  flight  owing  to  a  failure  in  the  write  circuitry  for  the 
tape  recorder.  Owing  Co  the  fact  that  the  automatic  data  reduction 
procedures  embodied  in  the  existing  computer  program  have  not  as  yet  been 
optimized,  these  data  are  considered  as  preliminary.  In  spite  of  this, 
however,  this  is  believed  to  represent  the  most  reliable  and  accurate 
source  of  information  on  stratospheric  particle  size  distribution 
presently  available. 

XI.  BRIEF  DESCRIPTION  OF  AEROSOL  COUNTER  SYSTEM 

The.  aerosol  counter  operates  on  the  principle  of  light  scattering, 
and  in  this  sense,  is  similar  to  other  instruments  of  this  kind.  The  novel 
feature  of  tnis  instrument,  however,  is  the  introduction  of  a  second  optical 
light  scattering  channel  which  compensates  ior  the  strong  dependence  of 
scattering  cross  section  on  index  of  refraction,  which  is  a  common  limita¬ 
tion  of  all  known  instruments.  Straightforward  solutions  of  the  Hie 
scattering  equations  has  shown  that  the  ratio  of  scattered  light  at  30 
and  10-  degrees  is  essentially  inrf<~,. «-ridcnt  of  the  index  of  refraction  over  the 
range  !.3  to  l . for  uen-wf  a'  ii<  spherical  particles.  M-'i'e'ver ,  this  ratio 
is  directly  proportional  to  the  particle  s lae  (expressed  in  the  forts  of 
the  scattering  parameter  a  •-  ed/>  ,  where  d  is  the  particle  diameter,  and 
X  is  the  wavelength  of  light  employed,  in  the  appropriate  units).  This 
greatly  relaxes  the  severity  of  assiuspt i ons  which  must  be  taada  in  infer¬ 
ring  particle  sise  i rocs  light  scattering  data.  Based  upon  this  result, 
the  device  was  designed  to  accoessodate  two  separate  annular  scattering 
channels,  viewing  the  scattered  light  in  the  S  to  15-degree  and  25  to 
35-degree  scat  ter ing  d irec t ions ,  respectively. 

Air  is  taken  into  the  chamber  by  a  specially-designed  sampling 
system  at  a  controlled  rate,  and  is  caused  to  flow  in  a  laminar  fashion 
through  a  small  bore  ( l  os  diameter)  tube,  which  is  broken  for  a  3  gee 
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portion  of  its  length.  Orthogonal  to  the  tube  axis  and  through  this 
broken  portion,  is  directed  an  intense  beam  (15nW)  of  highly  collimated 
light  at  632S&  from  a  Helium-neon  gas  laser.  The  collimation  feature 
is  most  important  since  it  allows  a  maximum  of  illumination  of  the  1  mm 
interaction  volume,  while  at  the  same  tune  maintaining  the  angle  of 
incidence  of  all  the  light  near  zero  degrees.  The  laser  is  maintained 
within  a  pressurized  housing  to  minimize  overheating  and  to  prevent 
high-voltage  arcing  at  the  reduced  ambient  pressures  encountered  during 
balloon  flight. 

Light  scattered  by  the  particles  within  the  interaction  volume  is 
collimated  by  two  annular  systems:  (1)  a  lens  for  the  10-degree  light, 
and  (2)  an  off-axis  parabaloid-of-revolution  mirror  for  the  30-aegree 
light.  The  two  collimated  beams  are  directed  to  each  of  two  PM  detectors 
by  a  set  of  two  45-degree  plane  mirrors  mounted  in  a  central  housing. 

The  detector  together  with  associated  condensing  lenses  and  image  aper¬ 
tures,  are  housed  in  two  lateral  arms.  The  image  apertures  are  carefully 
adjusted  to  allow  only  the  scattered  light  from  the  interaction  volume 
to  strike  the  PM  detectors.  Light  scattered  from  9uch  Boiid  surfaces 
as  the  end  of  the  flow  tubes  is  thereby  strongly  attenuated. 

Unused  light  from  the  laser  beam  passes  on  through  axially,  and  is 
absorbed  by  a  Wood's  horn  which  reduces  the  back  scattered  laser  light  in 
the  vicinity  of  the  interaction  volume  to  acceptable  levels.  Since  the 
system  sensitivity  is  limited  by  noise  associated  with  Rayleigh  scattered 
Light  by  ambient  air  molecules  within  the  small  interaction  volume,  all 
other  sources  of  laser  light  in  the  vicinity  of  this  region  must  be 
reduced  to  one  part  in  10^  which  is  accomplished  by  the  Wood's  horn 
arrangement 

The  flow  speed  is  maintained  constant  throughout  the  flight  by  means 
of  a  critical  orifice,  such  that  the  dwell  times  of  particles  crossing 
the  laser  beam  in  the  sensing  gap  is  about  250  microseconds.  Hence, 
pulses  of  this  duration  are  generated  at  each  of  the  two  PM  detectors, 
whose  peak  values  are  proportional  to  the  respective  scattered  light 
intensity  in  each  of  the  two  (10  and  30-degree)  scattering  channels.  The 
data  from  both  channels  are  digitized  and  recorded  on  a  Kinelogic  9 
channel  tape  recorder  in  a  format  compatible  with  a  IBM  360  computer. 

A  dynamic  range  of  four  orders  of  magnitude  is  accommodated  by  the  cir¬ 
cuitry  In  addition,  several  important  housekeeping  functions  are  also 
sampled  periodically,  and  recorded  on  tape  in  digitized  form.  Included 
here  are  PM  high  voltage,  laser  output  power,  air  flow  rate,  altitude, 
temperature,  and  other  parameters 

The  signal  processing  electronics  include  pulse  peak  detection 
circuitry,  clock  and  timing  circuits,  (which  generate  the  basic  timing 
waveforms  tor  the  entire  system),  read-in  control  circuits,  (which  control 
the  flow  of  information  traffic  through  the  sequence  of  complex  signal 
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processing  operations  being  performed  therein),  A/D  circuitry,  (which 
includes  a  unique  logarithmic/linear  function  to  handle  the  large  dynamic 
range  of  certain  signals),  and  magnetic  tape  recording  circuitry  (which 
provides  an  output  compatible  with  the  IBM  360  system,  including  parity, 

CRC  and  LRC  generators) . 

A  computer  program  has  been  generated  to  reduce  the  raw  data  and 
present  the  results  in  the  form  of  a  size  distribution,  together  with  the 
total  number  of  particles  being  counted. 

III.  PREPARATION  FOR  FIELD  MEASUREMENTS 

The  major  effort  under  the  current  contract  was  devoted  both  to  the 
repair  and  renovation  of  damaged  modules  and  the  fabrication  of  certain 
new  circuits  designed  under  the  previous  contract.  During  the  course  of 
this  work,  it  was  found  appropriate  to  modify  several  additional  modult s 
not  originally  contemplated.  As  part  of  this  activity,  the  equipment  was 
subjected  to  a  full  flight  simulation  test  at-  the  AVCO  Environmental  Test 
Facility  at  Wilmington,  Massachusetts. 

Upon  completion  of  these  renovations  and  modifications,  the  equipment 
was  calibrated  with  a  broader  range  of  particle  sizes  than  had  been  employed 
previously.  A  test  with  particles  having  a  different  index  of  refraction 
was  also  performed  which  demonstrated  that  the  equipment  is  much  more  in¬ 
sensitive  to  index  of  refraction  than  single-channel  devices. 

A-  RENOVATION,  MODIFICATIONS  AND  ENVIRONMENTAL  TESTING 

With  regard  to  renovation,  a  somewhat  unusual  procedure  was 
developed  as  a  result  of  the  rather  unique  circumstances  encountered  in 
the  present  situation.  Although  all  the  circuit  boards  had  been  dis¬ 
colored  (which  suggested  severe  deterioration  as  a  consequence  of  the 
test  chamber  accident  discussed  previously),  more  detailed  testing  of 
selected  circuits  indicated  that  operating  specifications  might  not  have 
been  compromised  to  us  serious  an  extent  as  had  originally  been  supposed. 

As  a  consequence,  it  was  decided  to  subject  all  circuit  boards  to  the 
wash,  rinse,  and  dry  cycle  of  an  automatic  washer  and  thereby  challenge 
each  element  severely.  A  detailed  microscopic  inspection  of  each  com¬ 
ponent  and  connection  was  then  performed ,  and  repairs  and/or  replacements 
were  executed  as  required.  This  was  followed  by  the  usual  testing  and 
checkout  procedure,  including  environmental  tests,  and  rework  if  necessary, 
etc.  This  procedure  was  predicated  on  the  judgment  that  the  resulting 
probability  of  component  failures  would  not  be  greater  chan  that  associ¬ 
ated  with  the  installation  of  new,  unchallenged,  components. 

The  following  subsections  describe  in  greater  detail  the  work 
performed  on  specific  items. 
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1.  Power  Supply  Renovation 

a.  Low  Voltage  and  Laser  Supplies 

These  two  modules  provide  power  for  operating  the 
electronic  processing  circuitry  and  laser  plasma  tube,  respectively,  and 
were  designed  and  operated  originally  as  two  separate  units.  The  low- 
voltage  supply  was  restored  to  operation,  but  an  adequate  assessment  of 
its  reliability  would  have  involved  a  greater  expense  than  that  required 
to  rebuild  the  unit  completely.  With  regard  to  the  laser  power  supply, 
a  means  was  found  to  greatly  simplify  its  design  and  hence  improve  its 
reliability.  In  the  course  of  implementing  these  modifications  it  was 
decided  to  combine  these  two  modules  into  a  single  dual  unit.  One  portion 
provided  a  22  kHz  400V  P-P,  square  wave  voltage  to  drive  the  laser  plasma 
tube.  This  was  current  limited  at  the  desired  load  level  to  control  the 
plasma  tube  current.  The  other  poition  provided  several  separate  outputs 
for  the  following  modules: 

6.3V  3.5  amps  -  pump  motor 

+  17V  -17V  -  electronic  circuitry 

+28V  (regulated)  -  tape  recorder  and  HV  power  supply 

b.  High  Voltage  Pqver  Supply  Renovation 

Preliminary  damage  inspection  and  testing  of  this 
module  indicated  that  renovation  could  be  accomplished  by  simply  replacing 
all  high  voltage  capacitors.  However,  when  renovation  was  initiated, 
several  more  subtle  aspects  of  damage  were  brought  into  clearer  focus. 
Among  these  was  deterioration  of  both  the  circuit  board  etched  wiring 
and  the  potting  compound,  as  well  as  the  possibility  of  latent  damage  in 
the  high  voltage  transformer  and  high  voltage  rectifier.  This  unit, 
pvocored  originally  from  a  separate  vendor  as  a  completely  encapsulated 
corafter icial  unit,  had  been  modified  previously  to  provide  short-circuit 
protection,  and  also  to  add  additional  filtering.  This  latter  feature 
had  been  incorporated  externally  to  avoid  the  extensive  repackaging  that 
would  have  been  required,  which  at  the.  time  was  precluded  by  the  existing 
scheduling  constraints.  Although  it  was  recognized  that  repackaging 
would  have  been  desirable  owing  to  the  geousetr  ical  constraints  peculiar 
to  this  particular  application,  it  was  never theioss  not  essential,  and 
therefore,  had  been  judged  to  be  of  secondary  importance  under  the  con¬ 
ditions  of  exigency  which  prevailed.  However,  such  was  not  the  case 
when  the  present  work  was  initiated.  The  need  to  perform  extensive 
repairs  rendered  it  a  relatively  simple  matter  to  repackage  the  high 
voltage  supply  with  the  added  filtering  and  the  short  circuit  protecting 
features  provided  as  an  integral  part  of  the  supply. 

The  high  voltage  power  supply  was  completed  and  tested 
thoroughly  in  the  AVCQ  Environmental  Test  Facility,  since  output  stability 
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was  an  important  consideration,  particularly  with  regard  to  temperature 
environment.  The  results  of  these  tests  showed  that  over  a  temperature 
span  of  -30°F  to  +150°F,  the  temperature  coefficient  of  the  output  voltage 
of  the  supply  was  well  within  acceptable  limits  averaging  approximately 
-0.015  percent/°F. 

2.  Tape  Recorder  Renovation 

The  tape  recorder  refurbishment  effort  included  the  replace¬ 
ment  of  four  of  the  mylar  drive  belts  and  the  replacement  of  the  inverter 
module.  Electrical  and  mechanical  tests  had  been  performed  on  the  tape 
recorder  at  an  earlier  date  to  establish  the  areas  requiring  rework.  Teats 
on  the  record  head  indicated  virtually  no  performance  degradation.  The 
motor,  on  the  other  hand,  showed  a  significant  increase  in  bearing  friction 
New  bearings  were  therefore  obtained  and  installed.  In  order  to  provide 
backup  support  for  the  field  operation,  however,  a  new  motor  was  procured 
as  part  of  the  spare  parts  inventory.  During  the  final  checkout,  it  was 
noticed  that  a  failure  had  occurred  in  one  of  the  buffer  amplifiers  which 
constitute  a  portion  of  the  input  circuitry  of  the  tape  recorder  unit. 

These  modules  were  found  to  be  redundant  since  the  output  commands  of  the 
aerosol  counter  electronics  had  been  designed  to  be  compatible  with  the 
tape  recorder  input  requirements.  Therefore,  rather  than  replace  the 
defective  amplifier,  it  was  decided  to  modify  the  tape  recorder  circuitry 
so  that  the  output  commands  would  be  applied  to  the  tape  recorder  inverter 
directly.  This  modification  eliminated  three  superfluous  buffer  ampli¬ 
fiers  and  hence  provided  improved  system  reliability. 

In  addition  to  the  components  discussed  above,  an  internal 
+5V  power  supply  is  incorporated  within  the  tape  recorder  module,  which 
constitutes  the  only  major  electronic  module  therein.  Since  this  supply 
showed  no  degradation  in  performance  and  is  very  lightly  loaded  with 
respect  to  its  output  current  capabilities,  it  was  concluded  that  a 
replacement  of  this  module  would  be  unwarranted. 

3 .  Ground  Support  Equipment 

A  rudimentary  ground  support  unit  had  been  fabricated 
previously  to  facilitate  testing  and  checkout  of  the  electronic  cir¬ 
cuitry  for  use  during  the  first  scheduled  field  exercises.  This  con¬ 
sisted  of  a  small  control  console  that  permitted  the  following  operations: 

(1)  use  of  either  internal  (i.e.  battery)  or  external 
power  for  the  entire  aerosol  counter 

(2)  selective  operation  of  the  tape  recorder  (on/off, 
fast /slow,  forward /reverse,  and  read/write) 

(3)  meter  monitoring  of  each  of  the  various  power 
supply  outputs 
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(4)  Monitoring  of  all  housekeeping  parameters  at  an 
accelerated  sampling  rate  to  permit  viewing  on  an 
oscilliscope,  or  at  the  normal  rate.  (During  normal 
operation,  these  parameters  are  sampled,  processed 
and  recorded  once  every  20  seconds ,  which  is  an 
inconvenient  sampling  rate  for  most  monitoring 
purposes . ) 

(5)  Activation  of  special  confidence  check  circuitry 
built  into  the  aerosol  counter.  (This  consists  of 

a  small  pulsed  light  source  to  simulate  the  scattered 
light  from  the  particles  observed  by  the  PM  tubes, 
which  constitutes  the  basic  datum  of  scientific 
import . ) 

(6)  Oscilliscope  monitoring  of  the  digital  data.  (This 
is  sampled  at  the  input  to  the  record  electronics 
module,  where  it  is  available  in  serial  form.)  Ap¬ 
propriate  signals  are  also  monitored  which  permit 
synchronization  of  the  oscilliscope  trace  to  facili¬ 
tate  viewing  of  the  digital  information. 

As  part  oi  _bn  damage  evaluation  and  preliminary  refurbish¬ 
ment  activity  under  the  prior  contract,  an  improved  ground  checkout  device 
was  designed  which  incorporated  all  of  the  above  features,  together  with 
a  digital- to-analogue  converter  and  associated  logic  circuitry  which 
permitted  the  monitoring  of  the  output  of  the  record  electronics  in  the 
form  in  which  the  information  is  impressed  on  the  tape  recorder  (i.e., 
parallel  format),  and  (2)  provision  for  monitoring  all  housekeeping 
functions  in  analog  form  before  digitalization.  These  combined  functions 
allow  a  comprehensive  check  to  be  performed  on  the  entire  digital  data 
processing  circuitry.  Additional  features  include  a  parity  check,  as 
well  as  provisions  for  feeding  out  several  other  synchronization  pulses 
to  facilitate  oscilliscope  monitoring. 

Under  the  present  program,  the  above-described  unit  was 
fabricated,  tested,  and  operated  in  the  field  exercises  with  entirely 
satisfactory  results. 

In  its  present  state  of  development,  the  ground  support 
unit  lacks  only  one  feature,  namely  the  ability  to  check  out  the  data 
following  impression  on  the  magnetic  tape.  Scheduling  and  funding 
limitations  precluded  the  incorporation  of  this  feature  under  this 
present  contract,  but  it  is  recommended  that  this  be  considered  for 
future  work. 


During  the  course  of  this  activity,  the  desirability  of 
having  g  battery  charger  which  would  be  provided  as  an  integral  part  of 
the  ground  support  equipment  was  recognized.  Although  the  batteries 
could  have  been  charged  externally,  this  would  have  resulted  in  the 
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inconvenience  of  either:  (1)  the  removal  of  the  battery  pack  to  the 
battery  charger  facility,  or  (2)  the  provision  of  a  separate  auxilliary 
power  supply  which  would  be  connected  to  the  ground  station  equipment 
by  an  external  cable.  In  order  to  provide  the  degree  of  operational 
flexibility  required  for  ground  checkout  during  a  field  operation,  a 
simple  battery  charger  circuit  was  built  and  incorporated  into  the 
system  as  in  integral  portion  of  the  ground  support  equipment. 

4,  Miscellaneous  Items 


a.  Elimination  of  Spurious  Pulses  in  CRC  Generator 

Detailed  visual  inspection  of  each  circuit  module  of 
the  refurbished  portions  of  the  aerosol  counter  was  performed  in  accordance 
with  the  program  plan.  After  the  environmental  test  had  been  performed 
at  the  AVCG  Test  Facility,  it  was  discovered  that  the  CRC  generator  had 
been  operating  intermittently.  This  was  traced  to  electrical  ringing  at 
the  trailing  edge  of  the  shift  pulse  for  the  CRC  circuitry,  and  was 
eliminated  by  adding  appropriate  additional  filtering. 

b.  Improvements  in  Confidence  Test  Circuitry 

As  a  result  of  the  use  of  an  AC  drive  signal  in  place 
of  the  original  DC  power  input  to  the  laser  exciter,  the  requisite  +150V 
source  of  power  was  no  longer  available  for  the  neon  lamp  light  sources 
in  confidence  test  circuitry  for  the  aerosol  counter.  Rather  than  design 
a  special  power  supply  for  these  lamps,  the  confidence  test  circuitry  was 
modified  to  work  with  light  emitting  diodes  (LED),  which  operate  with  low 
voltages.  These  have  been  available  commercially  only  recentl)  and, 
unlike  the  neon  bulbs  previously  used,  do  not  exhibit  a  "dark  effect". 

This  hat.  resulted  in  improved  system  reliability. 

c.  Peak  Detection  Circuits 

The  two  photomultiplier  preamplifier  modules,  using 
the  new  peak  detection  circuitry  designed  under  the  previous  contract, 
were  constructed  and  tested.  This  design  supplants  the  older  concept, 
in  which  that  portion  of  each  PM  output  pulse  protruding  above  a  fixed 
threshold  level  (determined  by  the  noise  level)  was  integrated  as  a 
measure  of  the  scattered  light  received  in  each  channel.  This  scheme 
had  been  incorporated  in  the  original  design  to  suppress  inaccuracies 
arising  from  noise  fluctuations  superposed  upon  the  pulses.  However,  in 
practice  this  was  not  found  to  be  a  significant  source  of  error,  so  that 
the  associated  disadvantage  (namely  the  introduction  of  inaccuracies 
associated  with  the  time  delay  of  the  initiation  of  a  pulse  emerging 
out  of  the  threshold  level),  r.u  longer  warranted  the  use  of  this  scheme. 
This  difficulty,  which  becomes  more  serious  the  smaller  the  particle  being 
observed,  has  been  alleviated  in  the  present  peak  detection  design,  since 
its  output  is  much  less  affected  by  pulse  initiation  time  delays. 
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d.  Time  Code  Generator 

The  overall  system  as  previously  constituted  was  found 
to  be  somewhat  deficient  with  regard  to  the  association  of  particle  size 
date  with  bt»llocn  altitude  data.  Although  the  requisite  information  was 
available  on  the  tape  record,  it  was  found  to  be  in  a  form  which  was  in¬ 
convenient  for  ready  data  processing.  Hence  to  facilitate  the  data 
reduction  procedure,  a  time  code  generator  was  designed,  fabricated  and 
incorporated  into  the  electronic  circuitry.  This  module  generates  the 
elapsed  time  in  hours  and  minutes  from  an  appropriate  reference  time 
(e.3.,  ’aunch)  and  records  this  information.  Two  vacant  tape  character 
positions,  which  had  previously  been  made  available  on  a  contingency 
basis  for  recording  housekeeping  parameters  which  might  be  found  neces¬ 
sary,  were  employed  for  this  purpose. 

In  this  connexion,  it  should  be  noted  that  correlation 
of  balloon  height  pro^'le  uata  between  the  three  separate  sources  (namely, 
aerosol  counter  altimeter,  balloon  control  system  altimeter,  and  radar 
retims),  exhibits  need  for  improvement.  Since  these  discrepancies  can 
be  traced  to  timing  inaccuracies  and  the  establishment  of  the  initial 
reference  time  betuee..  the  separate  tuning  systems,  it  is  recommended 
that  a  suitable  crystal-controlled  oscillator  be  installed  in  the  aerosol 
counter  as  the  clock  reference  of  the  time  code  generator,  and  that  a 
more  positive  means  be  instituted  for  establishing  an  intial  reference 
time  for  each  flight.  Timing  precisions  to  within  a  few  seconds  is  con¬ 
sidered  to  be  sufficient  for  most  balloon  flight  experiments.  It  is  of 
course  also  necessary  to  insure  that  similar  precision  is  available  on 
the  other  systems. 

e.  Additional  Temperature  Probus 

The  possibility  of  overheating  during  a  daytime  balloon 
flight  was  considered.  As  a  precautionary  measure,  the  gondola  was  covered 
with  specially  treated  aluminum  panels*.  These  panels  are  coated  with  a 
material  having  increased  emissivicy  in  the  1R,  to  reduce  the  sun  load. 
Nevertheless,  sufficient  uncertainty  remairad  regarding  the  predicted 
temperature  history  of  such  critical  component',  s,  for  instance,  the 
laser  housing,  power  supplies,  etc.,  to  warrant  the  use  of  additional 
temperature  monitors  wherever  possible.  The  de’ucion  of  the  corona 
altimeter  created  a  vacancy  in  one  of  the  tape  character  positions  in 
the  recorded  format  of  housekeeping  data.  1  1  is  permitted  the  inclusion 
of  a  separate  temperature  sensor,  for  which  the  appropriate  encoding 
circuitry  was  constructed  and  installed. 


Finished  by  ALZAK  process,  Alcoa  Aluminum  Corporation. 
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f,  Housekeeping  Monitor  Outputs 


The  analog  housekeeping  circuitry  used  in  the  original 
system  were  contained  in  three  circuit  boards.  These  were  replaced  by 
three  new  circuit  boards  in  which  the  design  has  been  modified  to  provide 
external  analog  outputs  in  parallel  form.  Thus,  continuous  housekeeping 
data  are  available  for  monitoring  the  input  to  the  A/D  converter.  This 
is  an  important  adjunct  to  the  ground  support  checkout  system. 

g.  Air  Pump 

A  new  air  sampling  pump  was  obtained  and  installed  to 
replace  the  pump  iv.  the  original  system  as  a  precautionary  measure  against 
the  possibility  that  its  life  might  have  been  shortened  as  a  result  of 
the  exposure  to  high  temperature. 


Environmental  Testing 


As  part  of  the  general  refurbishing  effort,  extensive  tests 
were  performed  on  the  electronics  portion  of  the  aerosol  counter,  using 
the  small  GCA  test  chamber  Temperature  tests  with  equipment  operating 
covered  the  range  from  -30°F  to  +150  F.  The  equipment  was  exposed  to 
several  cycles  of  temperature  stress,  which  served  the  useful  purpose  of 
identifying  not  only  marginal  components,  but  also  several  minor  circuit 
weaknesses  which  might  otherwise  not  have  been  detected. 


Full-scale  environmental  tests  of  the  complete  aerosol 
counter  system  were  then  conducted  at  the  AVCO  Test  Facility  in  Wilmington, 
Massachusetts.  These  tests  revealed  two  circuit  design  deficiencies  re¬ 
quiring  relatively  minor  modifications:  an  incorrect  scale  factor  in  the 
altitude  housekeeping  circuitry,  and  ringing  in  the  CRC  generator.  Cor¬ 
rective  measures  to  the  former  involved  changing  a  resistor,  while  those 
pertaining  to  the  latter  involved  the  addition  of  a  suitable  filtering 
network,  which  has  already  been  discussed  in  Subsection  4a  above.  Appro¬ 
priate  bench  tests  were  conducted  upon  completion  of  these  modifications, 
since  these  were  judged  not  to  be  extensive  enough  to  warrant  repetition 
of  the  full-scale  environmental  tests. 


The  optics  subassembly  was  installed  in  the  balloon  gondola 
and  checked  out  completely.  The  alignment  of  the  laser  and  the  10  and  30 
degree  optics  was  tested  and  adjusted  for  optimum  performance.  The  laser 
and  high-voltage  photomultiplier  housings  were  leak  tested  and  found  to 
be  pressure  tight.  The  sampling  chamber  was  found  to  have  several  small 
leaks  which  were  repaired. 


The  photomultipliers  and  laser  were  tested  in  an  environ¬ 
mental  temperature  chamber  over  the  range  -15°F  to  +100°F  and  were  found 
to  perform  satisfactorily.  Sensitivity  changes  were  similar  for  the  two 


tubes  and  the  overall  change  was  only  ~  10  percent,  which  is  tolerable. 
The  photomultipliers  exhibited  less  than  5  percent  hysteresis  in  sensi¬ 
tivity,  and  this  exhibited  partial  recovery  after  an  overnight  period. 
It  is  felt  that  full  recovery  would  be  realized  over  a  correspondingly 
longer  time  period. 


Air  Sampling  Flow  Rate 


The  system  was  designed  to  accommodate  static  air  sampling 
rates  of  3  to  5  cnr/sec.  This  range  is  associated  with  several  equipment 
parameters,  the  dominant  one  being  the  dwell  time  of  a  particle  in  the 
airflow  sensing  gap  as  it  crosses  the  orthogonal  laser  beam  (~1  mm  diam). 
This  represented  a  compromise  in  the  complex  trade-off  considerations 
between  large  sampling  rates  (to  increase  the  counting  statistics  for  the 
expected  low  aerosol  concentrations  at  balloon  altitudes),  a  small  diam¬ 
eter  flow  tube  (to  decrease  the  probability  of  more  than  one  particle 
being  illuminated  by  the  laser  beam  simultaneously),  and  a  slow  flow 
speed  (to  reduce  the  signal  processing  electrical  bandpass  and  thereby 
improve  signal-to-noise  ratio  for  small  particle  detection).  On  this 
basis,  the  electrical  bandpass  was  chosen,  and  the  circuitry  designed 
accordingly  to  accept  pulses  of  about  1/4  millisecond  minimum  duration. 
The  flow  rate  was  set  at  the  lower  end  of  this  range  to  favor  somewhat 
the  sensitivity  of  the  instrument  to  small  particles  at  the  expense  of 
higher  confidence  statistics. 


At  this  point,  a  complicating  feature  of  the  sampling 
system  should  be  pointed  out.  Due  to  the  necessity  of  avoiding  any 
windows  in  the  path  of  the  laser  beam,  (which  would  result  in  unacceptably 
high  background  light  levels  owing  to  surface  scattering  effects)  the 
sensing  gap  (~1  mnr)  communicates  with  a  much  larger  volume  (~7  liters) , 
which  is  sealed  from  the  surrounding  air  to  ensure  that  all  the  air  taken 
in  through  the  inlet  tube  indeed  passes  through  the  sensing  gap.  This 
volume  is  comprised  of  the  sampling  chamber  as  well  as  the  adjacent  laser 
baffling  chamber,  of  which  the  latter  is  by  far  the  larger  contributor. 


As  the  balloon  rises  through  the  exponentially  decreasing 
atmospheric  density,  the  air  in  this  chamber  must  be  evacuated,  and  can 
only  exit  through  the  sensing  gap.  Since  the  total  volume  flow  rate  of 
air  exiting  from  the  system  is  controlled  by  a  critical  orifice,  and  is 
therefore  constant,  the  air  sampling  rate  will  be  diminished  by  the 
chamber  evacuation  rate. 


Selection  of  the  original  flow  parameters  was  predicated 
on  a  balloon  rise  rate  of  100  meters /minute,  and  an  evacuation  volume  of 
about  1  liter.  On  this  basis,  calculations  indicated  (Ref. 2)  an  anti¬ 
cipated  maximum  evacuation  flow  rate  of  0.23  cnP/sec  which  was  less  than 
10  percent  of  the  static  air  sampling  rate,  and  hence  easily  tolerable. 
However,  the  estimates  of  balloon  rise  rates  and  evacuation  volume  were 
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found  subsequently  to  be  optimistic  by  a  considerable  margin.  More 
realistic  flight  profiles  indicated  that  balloon  rise  rates  would  likely 
reach  200  raeters/minute  during  certain  portions  of  the  flight.  In  addi¬ 
tion,  the  laser  baffle  design  which  finally  emerged  was  much  more 
voluminous  than  originally  anticipated,  resulting  in  an  evacuation 
volume  of  7  liters.  The  combined  effect  of  these  two  parameters  was  to 
increase  the  peak  evacuation  flow  rate  to  3.2  cm^/sec. 

To  reduce  the  possibility  of  negative  sampling  flow  (i.e., 
net  outflow  through  the  intake  port),  the  static  air  sampling  rate  was 
increased  to  5  cnr/sec.  This  was  achieved  by  simply  increasing  the  diam¬ 
eter  of  the  critical  orifice,  (the  pump  having  a  sufficient  capacity  to 
maintain  the  necessary  pressure  drop  across  the  constriction  under  the 
increased  flow  rate  to  assure  Mach  1  conditions  therein).  In  addition, 
the  differential  pressure-flow  transducer  was  recalibrated  to  incorporate 
the  higher  total  flow  rates  expected  during  descent,  which  result  from 
filling  flow  being  added  to  the  static  air  sampling  rate. 

In  spite  of  this  increased  value,  the  actual  balloon  flight 
profile  was  such  that  negative  sampling  flows  were  experienced  on  several 
occasions  during  balloon  rise.  These  portions  of  the  data  were  simply 
ignored.  It  should  be  noted  in  this  connection  that  the  particle  counting 
statistics  for  the  descent  portion  of  the  balloon  flight  are  much  more 
favorable  than  those  for  the  rise  portion,  owing  to  the  increased  total 
sampling  flow  under  these  conditions.  This  is  somewhat  fortuitous,  in 
that  the  descent  portion  data  are  less  subject  to  the  possible  disturbing 
influence  of  balloon  contamination  and  hence  are  likely  to  be  more  reliable 
by  virtue  of  this  fact  alone. 

B.  CALIBRATION 

As  noted  previously,  the  calibration  activity  included  two 
separate  particle  generation  schemes.  The  main  technique  involved  the 
use  of  a  suitable  nebulizer  and  highly  dilute  suspensions  of  commercially 
available  monodisperse  polystyrene-latex  spheres.  This  technique  had 
been  employed  previously  using  two  particles  sizes:  0,557  micron  and 
0.796  micron  and  is  described  in  Reference  1.  In  the  present  instance, 
seven  different  sizes  were  procured  (Refs.  3,4)  and  utilized:  0.198, 

0.312,  0.365,  0.500,  0.557,  0.796,  and  1.099  micron  diameter.  These 
particles  were  introduced  sequentially  into  the  aerosol  counter,  and  all 
data  were  processed  initially  as  though  these  sizes  were  unknowns  using 
an  assumed  calibration  constant.  The  output  for  each  particles  size  was 
then  compared  with  the  known  values,  and  a  corrected  average  value 
calibration  constant  was  thereby  obtained.  This  was  reinserted  into 
tfe  computer  program  in  iterative  fashion  to  converge  to  the  final  value  . 
...e  resulting  size  distribution,  plots  obtained  from  these  final  computer 
runs  showed  sharply  defined  peaks,  each  within  +10  percent  of  the  respec- 
tive  values  of  particle  sizes  noted  above. 

As  will  be  discussed  in  the  next  subsection,  a  different  numerical  calibra¬ 
tion  constant  applies  for  particles  with  cc  <  2  (i.e.  diam  <  0.428),  since 
the  absolute  value  of  scattered  intensity(as  opposed  to  the  ratio),  is 
utilized  in  this  range  of  particle  sizes. 
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The  polystyrene-latex  particles  employed  in  this  calibration 
scheme  have  an  index  of  refraction  of  1.58,  whereas  the  aerosol  counter 
was  intended  for  use  in  the  stratosphere  where  the  particle  index  of 
refraction  is  more  likely  to  be  in  the  range  1.33  to  1.50.  It  should 
be  noted  that  the  principle  upon  which  the  instrument  design  was  predi¬ 
cated,  (and  upon  which  the  interpretation  of  results  is  based)  embodies 
the  assertion  that  the  output  is  not  affected  by  the  index  of  refraction 
over  this  range, which  rests  on  sound  theoretical  reasoning.  Nevertheless, 
a  test  of  this  hypothesis  was  judged  to  be  highly  desirable.  Hence,  an 
alternative  scheme  was  devised  to  provide  particles  with  a  different 
index  of  refraction  than  that  exhibited  by  the  easily  available  polystyrene- 
latex  spheres. 

For  this  purpose,  aerosols  were  generated  using  DOP  (dioctal 
phtalate) ,  a  liquid  of  low  vapor  pressure,  in  combination  with  a  condensa¬ 
tion  type  aerosol  generator  (Refs.  5,6).  The  limited  scope  of  this  particu¬ 
lar  effort  permitted  the  generation  of  only  one  size  of  particles  for 
insertion  in  the  aerosol  counter.  The  average  diameter  was  measured  with 
a  specially-constructed  "owl"  device,  a  recent  version  of  which  is  de¬ 
scribed  in  Reference  5.  The  principle  involved  in  this  measurement  is  the 
observation  of  the  scattering  angle  at  which  the  first  red  band  of  the 
Tyndall  spectrum  appears.  Kitani  (Ref.  7)  describes  this  method  in 
detail  and  shows  that  over  the  range  of  indices  of  refraction  between 
1.33  and  1.55,  the  relationship  between  particle  radius,  r  and  the 
scattering  angle,  6  (from  the  forward  direction)  at  which  the  first  red 
band  appears  is  given  by: 

log(9/10)  +  1.43  log(lOr)  *  1.43  . 


Liu  (Ref.  5)  modified  this  equation  and  obtained  a  more  explicit  form: 


D 

P 


100 


-0.  7 


where  Dp  is  the  particle  diameter.  DOP  particles  have  an  index  of 
refraction  of  1.49  (Ref.  5)  and  are  therefore  well  suited  for  this 
purpose . 


A  single  particle  size  run  of  27-minute  duration  was  conducted 
using  this  scheme.  The  undiluted  output  of  the  evaporation-condensation 
generator  was  made  to  flow  through  the  "owl"  for  size  measurement,  which 
was  found  to  be  typically  0.9  micron.  A  small  fraction  of  the  concentrated 
aerosol  stream  from  the  generator  was  diluted  by  several  orders  of  magni¬ 
tude  with  a  clear  air  stream  so  that  the  concentration  entering  the 
aerosol  counter  would  be  sufficiently  small  to  preclude  the  simultaneous 
illumination  of  more  than  one  particle  by  the  laser  beam. 
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The  digital  output  of  the  aerosol  counter  as  recorded  on  magnetic 
tape  was  processed  using  the  same  procedure  as  would  be  employed  for  the 
actual  ballocn  flight  data.  The  results  were  shown  to  be  in  good  agreement 
with  the  measurements  obtained  using  the  "owl"  apparatus.  Thus,  it  can 
be  asserted  with  reasonable  assurance  that  the  particle  sizing  accuracy 
of  the  aerosol  counter  is  within  +10  percent  over  the  1.33  to  1.50  index 
of  refraction  range. 

Finally,  it  should  be  noted  that,  in  the  case  of  the  polystyrene 
spheres,  the  smallest  particles  tested,  namely  0.198  micron,  were  pro¬ 
cessed  without  difficulty.  From  visual  observation  of  the  PM  preamplifier 
outputs  of  one  of  the  two  scattered  light  channels,  the  pulses  observed 
were  well  above  the  average  noise  level.  This  would  imply  that  the  system 
should  be  able  to  approach  a  0.1-micron  capability,  at  least  under  the 
controlled  conditions  prevalent  in  the  laboratory. 

While  admittedly  this  is  a  rather  poorly  defined  criterion  upon 
which  to  extrapolate  system  behavior,  it  nevertheless  represents  the  best 
information  presently  available.  This  assertion  seemingly  conflicts  with 
conclusions  drawn  from  the  actual  flight  data,  as  will  be  discussed  below. 
However,  it  should  be  noted  that:  (1)  these  data  represent  the  only 
field  experience  with  this  equipment  to  data,  (2)  the  operating  conditions 
probably  deviated  considerably  from  those  prevailing  in  the  laboratory, 
and  most  importantly,  (3)  certain  data  processing  difficulties  have  been 
revealed  by  this  first  experience,  which  still  remain  to  be  rectified. 

It  is  therefore  anticipated  that  some  improvement  in  system  size  sensi¬ 
tivity  can  be  achieved,  which  will  narrow  the  gap  between  these  seemingly 
conflicting  conclusions.  Specific  recommendations  for  modification  of  the 
computer  processing  procedures  in  pursuit  of  this  goal  are  discussed  in 
the  last  section  of  this  report. 

C.  COMPUTER  PROGRAM  MODIFICATIONS 


The  program  in  its  original  form  contained  an  averaging  procedure 
which  selected  a  prescribed  set  of  altitude  increments  from  the  housekeeping 
data  and  processed  these  to  derive  average  value  size  distributions  for 
each  altitude  increment.  Early  in  the  contract,  consideration  was  given 
to  the  effect  of  likely  deviations  in  balloon-rise  flight  profile  from 
that  upon  which  the  original  averaging  procedure  was  predicated.  For  a 
smooth  rise  profile,  automatic  computer  processing  can  be  achieved  by 
relatively  simple  means.  However,  in  the  presence  of  significant  depar¬ 
tures  from  the  smooth  rise  case,  automatic  processing  of  the  data  would 
require  the  introduction  of  elaborate  look-ahead  and  smoothing  procedures, 
involving  not  only  extensive  revision  of  the  computer  codes,  but  also  other 
undesirable  complexities  such  as  extra  data  tapes  for  temporary  data 
storage,  and  much  longer  processing  times,  etc. 
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Upon  consultation  with  the  project  technical  monitor,  it  was 
decided  to  modify  tbe  averaging  procedure  whereby  time  (or  its  equivalent, 
namely  record  number)  is  used  as  the  independent  variable.  This  has 
resulted  in  some  simplification  of  the  computer  program.  However,  it  is 
no  longer  completely  automatic,  in  that  the  averaging  time  intervals  must 
be  inserted  externally  after  examining  the  actual  flight  profile.  This 
involves  a  slight  inconvenience,  namely  that  the  processing  cannot 
proceed  until  this  input  information  has  been  made  available.  However, 
based  upon  limited  experience  so  far,  this  does  not  seem  to  constitute 
a  serious  problem. 

Two  additional  modifications  were  introduced  as  a  result  of 
the  early  experience  gained  upon  initial  processing  of  the  actual  flight 
data.  The  first  of  these  modifications  was  occasioned  by  the  rapid 
changes  in  balloon  ascent  rate  experienced  on  the  first  balloon  flight, 
which  is  reflected  in  unexpectedly  low  flow  rates  in  the  particle  sensing 
gap.  As  a  result,  the  particle  size  distribution  determination  is 
rendered  unreliable  owing  to  the  possibility  of  a  single  particle  being 
counted  several  times  in  successive  sampling  intervals.  In  addition, 
the  total  particle  count  in  these  instances  is  subject  to  large  statistical 
uncertainties  owing  to  the  small  number  of  counts  involved.  This  is  also 
aggravated  by  the  fact  that  the  total  time  duration  for  a  given  altitude 
increment  is  reduced  with  increasing  rise  rate.  For  these  reasons  the 
program  was  modified  to  ignore  all  records  with  small  (or  negative)  flow 
rates . 


The  second  of  these  modifications  was  concerned  with  that 
portion  of  the  program  which  yields  the  total  concentration.  Since  in 
this  case  the  larger  number  of  counts  involved  result  in  smaller 
statistical  uncertainties,  it  was  possible  to  modify  the  program  to 
enab’e  the  total  particle  concentration  to  be  extracted  from  a  smaller 
number  of  records.  This  permitted  these  particular  data  to  he  presented 
on  a  finer  altitude  grid. 

The  implementation  of  this  latter  modification  introduced  a 
minor  problem,  namely  that  record  counting  is  upset  if  the  final  record 
which  is  processed  in  this  particular  manner  for  a  given  altitude  incre¬ 
ment,  happens  to  be  a  record  with  low  or  negative  flow  rate.  This 
problem  is  a  direct  consequence  of  the  fact  that  these  latter  records 
are  now  ignored  by  virtue  of  the  previous  modification  discussed  above. 
This  was  not  recognized  as  such  until  after  the  data  processing  had  been 
completed.  Since  it  was  not  possible  to  make  the  necessary  adjustments 
prior  to  the  completion  of  work  under  the  present  contract,  this  problem 
remains  to  be  rectified. 
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IV  FLIGHT  TESTS 

A.  FIELD  EXERCISES 


The  field  operation  is  discussed  in  somewhat  greater  detail 
than  would  ordinarily  be  the  case  for  a  report  of  this  kind,  owing  to 
the  fact  that  this  represents  the  initial  field  experience  with  the 
system. 


The  instrument  was  shipped  by  air  to  Holloman  Air  Force  Base, 
and  checked  out  for  a  second  series  of  environmental  tests,  which 
included  an  altitude  run  (i.e.  without  temperature  variation)  on 
15  October  1970,  and  a  complete  alt itude/teraperature  simulation  profile 
on  16  October.  During  the  course  of  the  latter  test,  for  reasons  beyor.d 
our  control,  the  instrument  was  placed  accidently  in  direct  contact  with 
the  refrigeration  coils  of  the  chamber,  so  that  the  low  temperature 
exposure  was  far  more  severe  than  anticipated  in  the  system  design.  In 
spite  of  this,  •Tie  equipment  performed  satisfactorily,  with  the  exception 
of  one  housekeeping  unit,  namely  the  elapsed  time  indicator  which  was 
repaired  subsequently. 

On  19  October  the  equipment  was  shipped  to  the  Stallion  Site 
of  the  White  Sands  Missile  Range  for  prelaunch  readiness  operations, 
which  were  performed  during  the  ensuing  week.  One  major  problem  area 
remained  to  be  investigated,  namely,  the  possibility  of  thermal  danger 
to  the  equipment  resulting  from  the  simultaneous  effects  of  internal 
power  dissipation  and  external  sun  load  at  high  altitudes,  where  heat 
transfer  by  conduction  to  the  atmosphere  would  be  diminished  signifi¬ 
cantly  by  virtue  of  the  reduced  ambient  pressure.  Owing  to  the  short 
time  available,  the  investigation  was  limited  to  a  determination  of  the 
relative  absorptivity  of  various  e  at ings  and  materials  that  could 
conceivably  be  used  to  replace  the  spec ially-coated  aluminum  sheets 
discussed  previously. 

Several  different  samples  were  exposed  to  direct  sunlight  on 
a  relatively  calm  day.  The  temperature  rise  of  each  sample  was  recorded 
after  the  teo^perature  had  attained  a  stable  value.  The  samples  included 
a  panel  of  the  type  used  to  cover  the  aerosol  counter  gondola,  an  aluminum 
panel  which  was  painted  with  a  white  titanium  dioxide  paint,  and  an  alumi¬ 
num  panel  which  had  been  iridited.  The  temperature  rise  of  the  gondola 
panel  was  measured  for  both  the  shiny  surface  and  the  dull  surface.  The 
test  results  showed  that  the  lowest  temperature  rise  was  experienced  by 
the  gondola  panel  for  the  shiny  side.  The  next  lowest  temperature  rise 
corresponded  to  the  gondola  panel  on  the  dull  side.  The  temperature 
rise  of  the  white  panel  was  significantly  greater  than  the  temperature 
rise  of  the  gondola  panel  even  when  compared  to  the  dull  side  The 
temperature  rise  of  the  iridited  aluminum  panel  was  the  greatest. 
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Because  of  both  the  relatively  crude  nature  of  the  experiment, 
and  other  factors  such  as  variations  caused  by  occasional  gusts  of  wind 
and  the  lack  of  appropriate  thermal  measuring  instrumentation,  it  was 
nui.  possible  to  quantitatively  determine  the  absorptivities  and  the 
emissivities  of  the  samples  to  a  useful  degree  of  accuracy.  Nevertheless, 
the  tests  served  to  corroborate  that  the  panels  which  had  already  been 
installed  on  the  gondola  were  the  most  appropriate  for  minimizing  tempera¬ 
ture  rise  therein.  Although  it  was  indicated  that  slightly  better  results 
would  have  been  obtained  had  the  panels  been  installed  with  the  shiny 
side  facing  outward,  it  was  felt  that  tha  difficulties  involved  in 
inverting  the  panels  would  not  be  worth  the  relatively  small  improvement 
to  be  attained  thereby. 

During  the  course  of  this  investigation,  an  opportunity  was 
presented  for  the  first  time  to  operate  the  entire  system  continuously 
for  several  hours  under  full  solar  illumination.  At  the  end  of  this 
period,  it  was  noted  that  laser  power  h3d  decreased  significantly  and 
that  the  laser  temperature  had  reached  65°C.  It  is  believed  that  pre¬ 
viously,  this  temperature  had  never  exceeded  55°C.  Inspection  of  the 
laser  and  subsequent  tests  indicated  that  the  laser  power  loss  was  due 
to  optical  misalignments  which  appeared  to  be  a  function  of  both  temper¬ 
ature  and  operating  time.  Further  investigation  disclosed  that  the 
mirrors  of  the  laser  are  joined  to  the  laser  tube  by  means  of  metal 
bellows.  These  bellows  act  as  a  shield  against  dust  deposits  on  the 
mirrors  and  Erester  windows  of  the  laser  tube,  and  hence,  are  r. latively 
air  tight.  The  laser  as  packaged  in  the  aerosol  counter  is  sealed  in  a 
tube  which  is  maintained  internally  at  local  atmospheric  pressure.  As 
the  laser  temperature  increases,  the  air  volume  inside  the  tube  is  heated 
and  the  pressure  therein  is  increased  enough  that,  by  virtue  of  the 
bellows  arrangement  and  the  nature  of  the  mirror  mount  design,  a  mis¬ 
alignment  develops.  To  correct  this  problem,  the  laser  housing  was 
fitted  with  a  pressure  value  and  gauge  to  allow  selection  of  an  optimal 
operating  pressure  just  prior  to  flight.  The  gauge  also  provides  an 
extra  assurance  that  Che  seal  of  the  laser  tube  has  no  leaks.  Such  a 
gauge  would  be  desirable  in  tne  high  voltage  section  which  also  is 
maintained  at  atmospheric  pressure. 

On  the  morning  of  28  October  1970,  the  aerosol  counter  was 
ready  for  its  initial  balloon  flight  test.  At  7:02  AN  the  package  was 
launched  from  the  runway  at  the  Stallion  Site  of  White  Sands  Missile 
Range.  The  progransoed  flight  profile  consisted  of  an  ascent  phase  which 
lasted  2-3/4  hours.  This  was  followed  by  a  float  phase  at  an  altitude 
of  about  80,000  feet  for  a  time  interval  of  about  20  minutes.  During 
the  final  three  hours  of  flight,  the  balloon  descended  slowly  to  about 
50,000  feet  at  which  point  a  telemetered  ground  command  caused  the 
balloon  to  be  separated  from  the  aerosol  counter  and  destroyed.  The 
experimental  package  was  then  returned  to  earth  by  parachute. 
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Unfortunately,  the  flight  lasted  somewhat  longer  than  the  automatic  turn 
off  time  of  6  hours,  so  that  no  data  were  col?^cted  during  the  latter 
portion  of  the  descent  phase  (below  about  16  km). 

The  gondola  landed  in  flat  grazing  terrain,  in  an  upright  posi¬ 
tion.  Impact  was  nearly  vertical  although  the  crush  pads  indicated  some 
horizontal  component.  A  limited  on-site  inspection  indicated  that  the 
equipment  had  probably  operated  as  planned,  and  data  had  been  obtained. 

The  tape  recorder  was  removed  and  a  dump  of  the  data  was  obtained  at  the 
White  Sands  Missile  Range  on  the  following  day.  It  was  intended  ori¬ 
ginally  that  the  data  would  also  be  processed  at  White  Sands;  however, 
the  center  computer  facility  personnel  had  not  had  time  to  program  the 
necessary  input  and  format  statements.  A  limited  manual  data  reduction 
was  performed  and  the  following  information  was  obtained: 

(1)  The  balloon  trajectory  contained  many  points  at  which 
rapid  changes  in  rate  of  rise  occurred.  Furthermore,  the 
rate  of  rise  at  many  points  was  considerably  higher  than 
the  scheduled  600  feet  per  minute  value.  Because  of  these 
two  undesirable  factors,  the  fluctuations  in  flow  rate  into 
the  sampling  volume  were  often  quite  severe.  In  fact,  at 
several  points  just  before  float  altitude  was  reached  the 
input  sampling  rate  dropped  to  zero. 

(2)  A  count  was  made  of  the  number  of  particles  which  were 
present  in  several  representative  records  on  the  magnetic 
tape  record.  A  crude  computation  was  then  performed  to 
obtain  a  particle  density  figure  at  several  altitudes  in 
units  of  particles  per  cubic  centimeter.  The  values  thus 
obtained  appeared  to  agree  with  the  currently  accepted 
estimates . 

(3)  Housekeeping  data  appeared  to  be  normal  except  for  the 
elapsed  time  output. 

(/•)  Laser  power  remained  relatively  stable  during  the  flight. 

(5)  A  relatively  large  number  of  parity  errors  were  indicated 
by  the  computer  printout.  These  could  result  either  from 
excessive  dust  on  the  magnetic  tape  or  in  an  intermittent 
noise  susceptibility  in  the  equipment.  It  is  also  quite 
possible  that  noise  introduced  by  the  power  supply  inverter 
could  have  caused  some  of  the  errors. 

(6)  On  at  least  two  files  of  data  the  computer  printout  indi¬ 
cated  that  no  end -of -file  mark  was  found  in  the  recording. 

The  equipment  difficulties  listed  in  Items  3  and  6  above  are 
relatively  minor  and  have  virtually  no  effect  on  the  quality  of  the  data 
collected.  Item  5  on  the  other  hand  could  result  in  the  loss  of  all  data 
in  several  records.  However,  this  is  not  a  particularly  serious  circum¬ 
stance,  since  only  about  30  percent  of  the  record*,  are  thus  affected,  and 
of  these  only  a  small  fraction  will  be  discarded  in  fact. 
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In  this  connection  it  should  be  noted  that  the  total  number  of 
successful  sampling  intervals  recorded  during  the  six  hour  flight,  from 
which  useful  data  can  be  extracted,  exceed  10  .  Thus,  despite  these  few 
minor  equipment  and  operational  problems  noted  above,  the  first  balloon 
flight  of  28  October  was  highly  successful.  The  data  have  been  processed 
in  a  preliminary  fashion  and  the  results  are  presented  in  Subsection  C 
be  low . 


On  the  morning  of  30  October,  the  aerosol  counter  was  returned 
to  the  Stallion  Site  and  checkout  preparation  for  the  second  balloon 
flight  was  conmenced.  The  equipment  was  found  to  be  operating  within 
specifications  except  for  (1)  the  end-of-file  circuit  and  the  readout 
portion  of  the  housekeeping  circuit  which  is  used  to  generate  elapsed 
time,  (2)  spurious  pulses  at  the  frequency  of  the  power  supply  inverter 
which  appeared  on  several  command  leads  to  the  tape  recorder  electronics 
circuitry,  and  (3)  the  output  voltage  of  the  5V  supply  for  the  integrated 
circuitry,  which  was  below  the  minimum  acceptable  level.  The  spurious 
pulses  were  eliminated  by  the  addition  of  appropriate  by-pass  capacitors, 
and  che  output  of  the  5V  supply  was  adjusted  to  the  appropriate  level. 
The  combination  of  these  two  corrective  measures  resulted  in  the  elimina¬ 
tion  of  the  end-of-file  circuit  difficulty. 


During  the  remaining  days  prior  to  launch,  several  additional 
equipment  difficulties  were  brought  to  light,  among  which  were  included: 
a  cracked  battery  cell,  laser  power  supply  transformer  failure,  and  a 
leak  in  t he  high  voltage  pressurized  compartment .  The  necessary  repairs 
were  all  executed,  and  the  entire  system  was  readied  for  launch  at  7:30  AM 
on  3  November  1970. 


As  in  the  case  of  the  first  flight,  the  experiment  package  was 
recovered  with  no  apparent  damage.  The  position  of  tape  in  the  tape 
recorder  indicated  that  the  pre-programmed  6-hour  operating  schedule  had 
taken  place.  However,  when  a  caaputer  playback  was  attempted  on  che 
following  day,  it  was  discovered  that  the  portion  of  the  tape  which 
corresponded  to  the  flight  was  blank.  The  only  data  on  the  tape  corre¬ 
sponded  to  the  initial  preflight  test  activity,  which  had  beer  obtained 
in  the  balloon  building  at  about  3:00  AM  on  the  day  of  che  flight.  Hence, 
the  failure  must  have  occurred  between  final  checkout  and  the  commencement 
of  actual  in-flight  re-ording. 

The  cause  of  the  failure  to  record  has  not  yet  been  identified 
positively.  The  roost  probable  cause  appears  to  be  the  loss  in  dc  power 
in  the  record  portion  of  the  tape  recorder  electronics  circuitry.  Sub¬ 
sequent  tests  indicate  that  the  system  is  recording  properly.  This 
points  to  an  interroittent  problem  which  must  be  investigated  and  cor¬ 
rected  prior  to  future  field  test  applications. 


19 


!  >4W 


WP3. 


B.  DATA  PROCESSING 
1.  Introduction 


As  noted  previously,  it  has  not  been  possible  to  develop 
optimum  data  reduction  procedures  within  the  limited  tine  and  funding 
available  under  the  present  contract.  Thus,  the  results  to  be  presented 
in  this  report  are  considered  as  preliminary,  pending  the  outcome  of 
several  experimental  computer  runs  which  remain  to  be  performed.  It  is 
anticipated  that  these  data  processing  improvements  will  be  reflected  in 
the  realization  of  the  full  potential  of  this  instrument,  namely  the 
attainment  of  0.1  to  0.2  micron  diameter  capability. 

A  description  of  the  manner  in  which  the  information  is 
sampled  and  recorded  on  tape  is  given  first,  together  with  pertinent 
comments  and  several  necessary  input  relationships  for  the  data  reduction 
computer  program. 

Also  discussed  is  the  method  by  which  background  levels 
are  sampled  and  pr.es  ted,  as  well  as  low  laser  power,  "dead"  time  cor¬ 
rections,  and  tempetacure  corrections,  are  to  be  applied. 

uis  is  followed  by  a  section  describing  the  rationale 
upon  which  the  automatic  computer  processing  of  the  scientific  data  is 
based. 


2.  Data  Reduction  Input  Considerations 
a.  Tape  Format 

The  recordings  prepared  by  the  aerosol  counter  cir¬ 
cuitry  consist  of  a  series  of  data  groups  called  records,  which  are 
separated  by  blank  areas  of  approximately  0.6  inch  in  length.  A  pic¬ 
torial  representation  of  a  typical  data  record  is  shown  in  Figure  1. 

Eight  separate  tracks  (labelled  0-7  but  not  in  numerical  sequence)  are 
used  for  recording  data,  while  the  ninth  track  (labelled  P)  is  used  for  a 
parity  bit.  As  can  be  seen  from  the  figure,  the  recording  begins  with  a 
16  character  unique  code  followed  by  16  characters  representing  house¬ 
keeping  functions.  At  the  end  of  the  first  32  characters  the  recording 
system  automatically  switches  from  a  housekeeping  mode  to  a  data  mode. 

Thus  the  remaining  3488  characters  are  employed  for  data  recording. 

Since  two  characters  are  required  for  each  data  sample  (10°  PM  output 
and  30°PM  output),  a  total  of  1744  samples  can  be  stored  in  each  record. 

(1)  The  first  eight  characters,  which  are  recorded  in 
positions,  (or  bytes)  0  through  7  of  a  record  con¬ 
tain  the  binary  representation  of  255  (all  bits 
are  "11 's").  The  parity  bit  is  a  one  because  odd 
parity  is  used. 
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Figure  L.  Pictorial  representation  cf  typical  data  record. 


(2)  The  following  eight  characters  (positions  8  through 
15  on  the  tape)  represent  0  in  binary  form  so  that 
the  parity  bit  for  these  characters  is  a  one. 

(3)  The  17th  character  (position  16)  through  the  30th 
character  (position  29)  contain  housekeeping  param¬ 
eters  in  straight  binary  form.  The  system  normaliza¬ 
tion  is  such  that  all  zeros  represent  zero  volts 

and  all  "l's"  represent: 

5V  x  (1  -  —t)  a  4.980  volts 
Thus,  the  voltage  representation  is 

EH  =  256  X 

where  Ag  is  the  binary  character  stored  on  the  tape. 
Tha  voltage  representation  Ejj  for  each  housekeeping 
parameter  can  be  translated  to  the  actual  parameter 
using  the  scale  factors  and  calibration  referenced 
in  Table  1. 

(4)  The  characters  in  positions  30  and  31  represent 
binary  time  in  hours  and  minutes  respectively. 

(5)  The  characters  which  are  present  in  positions  32 
through  position  3519  represent  samples  recorded 
during  the  data  mode  interval.  The  32nd  and  33rd 
positions  contain  the  first  data  sample.  The  even 
numbered  positions  from  this  point  represent  the 
output  of  the  photomultiplier  at  the  10  degree 
scattering  angle  while  the  odd  numbered  positions 
represent  the  PM  output  at  30  degrees.  In  this 
sample  and  in  subsequent  data  samples  the  character 
positions  are  such  that  if  N  represents  the  char¬ 
acter  position  on  tape  of  the  10  degree  sample,  the 
position  of  the  corresponding  30  degree  channel  will 
be  at  N  +  1  (never  at  N  -  1).  During  the  data  mode, 
the  binary  number  in  each  character  position  is 
related  to  the  output  of  its  corresponding  photo¬ 
multiplier  output  as  follows: 

-S  2  Ar  +  1 

Ip  =  10  x  10  exp  (  £28  ^  x  (l  +  0.0183) 
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where  Ip  is  the  associated  photomultiplier  current 
and  Ab  is  the  binary  character  stored  on  the  tape. 
For  Ab  =  255,  Ip  =  0.S82  x  10~^<1  +  0.018)  which 
corresponds  to  the  maximum  anticipated  level  of 
input  current  from  the  photomultiplier  anode.  The 
minimum  photomultiplier  output  occurs  at  Ag  =  0 
which  corresponds  to  Ip  <  1.037  x  10~8  amperes. 

At  the  end  of  the  3520th  character,  the  special 
GRC  character  is  inserted  after  a  four  character 
blank  interval.  Finally,  the  LRC  character  is 
recorded  four  additional  characters  later  followed 
by  a  0.6  inch  blank  interval  representing  the  inter¬ 
record  gap. 

b.  Sampling 

The  correction  factors  required  as  a  result  of  the 
method  of  sampling  are  now  discussed.  The  method  of  sampling  is  such  that 
time  is  divided  into  independent  increments  of  0.01  second.  During  the 
first  11/16  portion  of  each  time  increment,  a  sample  gate  is  in  the  active 
state  and  the  output  of  the  PM  tubes  is  monitored  for  the  detection  of  an 
output  signal.  The  sampling  interval  is  followed  by  the  readout  interval 
which  takes  place  during  the  remaining  5/16  portion  of  tne  0. 01-second 
time  increment.  The  readout  process  consists  of  converting  the  outputs 
of  the  two  photomultiplier  channels  to  digital  form  and  translating  the 
digital  output  to  a  form  suitable  for  preparing  an  IBM  compatible 
recording. 


c.  Method  of  Correcting  for  Background  Level 

If  no  signal  pulses  are  detected  in  the  11/16 
sampling  interval  of  time  increment,  a  250-microsecond  time  strobe  is 
generated  beginning  250  microseconds  before  the  end  of  the  sampling 
interval.  This  permits  a  sample  to  be  taken  of  the  background  output 
level  of  the  two  PM  tubes  for  every  time  increment  during  which  no  signals 
are  detected.  Thus  an  output  will  be  available  for  processing  at  the  end 
of  every  sampling  interval  whether  an  aerosol  particle  is  detected  or  not. 
The  recording  of  background  levels  is  useful  because  it  allows  the  signal 
contributions  resulting  from  Rayleigh  scattering,  dark  current,  dc  shifts 
in  the  electronics  circuitry,  etc.  to  be  subtracted  out  during  the  data 
reduction  process. 


d.  Correction  for  Sampling  "Dead"  Time 

Whenever  a  particle  is  detected,  the  sampling 
process  is  turned  off  for  the  balance  of  that  0.01  second  time  increment. 
Thus,  only  the  first  particle  which  is  detected  in  each  interval  will  be 
measured  and  recorded.  In  the  computer  programming,  a  correction  must  be 
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made  both  for  the  5/16  inactive  interval  in  each  time  increment  and  for 
the  inactive  interval  following  the  detection  of  each  particle.  It  should 
be  noted  that  unlike  the  Geiger  counter  which  has  a  fixed  "dead"  time 
following  each  detected  event,  this  equipment  has  a  variable  "dead"  time 
depending  on  which  portion  of  the  0.01-second  time  increment  an  event 
occurs.  The  multiplying  factor  by  which  the  data  rates  are  corrected  is 
given  by 

16  ^n(Nfc/N  ) 

p  =  Ti  (i-N  /n  5" 

where  N0  is  the  number  of  intervals  out  of  Nt  successive  intervals  during 
which  no  pulses  are  detected. 

Because  of  the  background  sampling  which  occurs  in 
every  sampling  interval  during  which  no  pulse  is  observed,  the  quantity 
N0  is  not  obtained  readily.  It  would  have  been  desirable  to  indicate  the 
occurrence  of  a  background  sample  by  means  of  a  flag  bit  associated  with 
each  character  recorded  on  tape  but  unfortunately  the  entire  8  bits  had 
to  be  used  to  meet  the  system  accuracy  requirements.  Thus,  it  was  not 
possible  to  incorporate  such  a  flag  bit  in  the  present  design.  Conse¬ 
quently,  it  is  necessary  to  determine  NQ  using  a  computer  program  to  scan 
the  data.  The  programming  should  be  such  that  only  data  which  corre¬ 
spond  to  light  pulses  with  amplitudes  less  than  the  system  triggering 
level  will  be  entered  into  the  N0  count. 

The  maximum  length  of  time  over  which  the  count,  Nt, 
may  be  taken  is  related  to  how  rapidly  the  statistical  properties  of  the 
particles  being  measured  are  changing.  If  the  count  rate  of  the  detected 
particles  is  not  changing  the  interval,  Nt,  may  be  set  arbitrarily.  For 
convenience,  the  number  of  samples  in  a  record  (Nt  =  1744)  might  be  used 
for  the  determination  of  the  ratio  N0/Nt. 

The  sampled  volume  associated  with  each  record  is  given 
by  1744  samples  x  0.01  second  per  sample  x  flow  rate  =  17.44  x  flow  rate. 

The  first  0.16  second  and  the  last  2.4  seconds  of  each 
20-second  record  interval  should  be  disregarded  because  these  intervals 
are  employed  for  housekeeping  and  the  inter-record  blank  interval  on  the 
tape.  During  these  two  intervals  any  particles  which  happen  to  be  de¬ 
tected  are  discarded  automatically. 

e.  Other  Computational  Considerations 

The  corrections  for  variations  in  volumetric  sampling 
rate  can  be  obtained  conveniently  by  making  use  of  the  flow  rate  parameter 
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which  is  recorded  during  the  housekeeping  mode  of  operation  when  these 
15  separate  parameters  are  sampled  and  recorded.  The  sampling  rate,  once 
every  twenty  seconds,  is  comfortably  higher  then  the  anticipated  rate  of 
change  of  any  of  the  housekeeping  parameters  during  the  experiment. 

An  absolute  value  of  laser  power  is  required  for  the 
smaller  particle  sizes  (below  a  =  2)  for  which  the  computation  makes  use 
of  the  direct  outputs  of  the  10  and  30-degree  channels  rather  than  the 
ratio  of  the  outputs.  It  is  anticipated  that  the  laser  power  parameter 
can  be  utilized  directly  in  the  computer  program  without  modification, 
'lowever,  because  the  laser  power  parameter  is  monitored  by  a  silicon 
photodiode  a  temperature  correction  may  be  required.  This  can  be  accom¬ 
plished  using  the  laser  temperature  parameter  which  is  available  as  one 
of  the  housekeeping  parameters  and  can  therefore  be  incorporated,  if 
needed,  into  the  computer  program.  The  temperature  parameter  is  non¬ 
linear;  therefore,  the  calibration  curve  which  is  provided  with  each 
temperature  sensor  is  used  in  the  determination  of  temperature. 

In  most  cases  the  accuracy  to  which  temperature  must 
be  known  in  order  to  perform  corrections  for  temperature  dependent  param¬ 
eters  is  not  very  critical.  Therefore,  in  place  of  a  calibration  curve  an 
equation  with  three  calibration  constants  (C^,  C2  and  C3)  can  be  used. 

In  terms  of  the  binary  character,  Ag>  stored  on  the  tape  the  equation  is 
given  by 


T  = - 

1  +  C2  ^n 


256  -  h  '  aBo  \ 

256  'ho  \  ) 


C 
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The  binary  character  Ag0  represents  the  binary  character  recorded  at  a 
temperature  of  25°G .  Normally  AgQ/ (256-Ag0)  %  1. 

Using  typical  values  of  the  calibration  constants  the 

equation  becomes 


T  = 


298 


,1  +  0.0806  2n 


(256-Ag) 

■(256^) 
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1  +  0.0806  jjn 


273°C 


Except  for  the  altitude  parameter,  the  remaining 
housekeeping  parameters  are  for  diagnostic  purposes  and  do  not  enter  into 
the  computations.  However,  for  the  first  few  balloon  flights,  the  com¬ 
puter  should  be  programmed  to  read  out  the  values  of  these  parameters 
(at  least  on  a  once-per-ten-minute  basis). 
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f.  IBM  Compatibility 


The  recording  format  is  such  that  compatibility  to  the 
IBM  360  standards  is  maintained  except  for  work  length.  In  this  system, 
the  eight  data  bits  of  each  character  recorded  on  tape  represent  an 
independent  eight  bit  binary  number.  The  ninth  bit  represents  the  ver¬ 
tical  parity  bit  and  will  be  a  "one"  whenever  the  other  eight  bits  con¬ 
tain  an  even  number  of  ones . 

3 .  Rationale  for  Processing  Scientific  Data 

The  primary  scientific  data  input  consists  of  pairs  of 
maximum  intensities  measured  in  the  10  and  30-degree  channels  for  a 
specific  sampling  time.  These  data  are  compared  with  those  derived 
from  Mie  scattering  calculations  as  modified  by  experimentally  deter¬ 
mined  calibration  data.  Using  these  calibration  data,  two  comparison 
curves  are  constructed,  the  first  corresponding  to  analysis  based  on  the 
ratio  of  the  scattering  intensities  into  the  10  and  30-degree  channel, 
which  is  appropriate  for  particles  with  a  exceeding  2.0  as  shown  in 
Figure  2.  The  second  curve  presented  in  Figure  3  corresponds  to  analysis 
based  on  the  absolute  intensity  of  the  scattered  signal  into  the  10-degree 
channel,  which  is  appropriate  for  particles  with  a  less  than  2.0. 

Care  must  be  exercised  to  obtain  continuity  at  the  (2-value 
of  2.0  This  is  a  natural  consequence  of  the  fact  that  two  different  data 
reduction  schemes  are  employed  at  higher  and  lowera -values  respectively, 
so  that  some  residual  discontinuities  arising  from  small  errors  in  the 
separate  calibration  constants  and/or  extrapolation  of  the  fitted  curves, 
are  to  be  anticipated. 

The  first  step  in  the  analysis  of  a  data  pair  consists 
of  making  a  decision  as  to  whether  it  is  likely  that  a  particle  was 
detected  during  the  sampling  period.  Let  I^q  and  I3C  rapresent  the 
recorded  signals  and  let  T^q  and  T30  represent  threshold  criteria  speci¬ 
fied  as  input  data.  One  has  four  possible  conditions: 

(1)  lj.0  <,  T10  and  I3Q  ^  T3Q.  This  condition  indicates 
that  no  particle  has  been  detected.  The  values  I^q 
and  I3Q  are  used,  along  with  other  pairs  in  the  record 
which  satisfy  the  same  condition  to  calculate  a  mean 
background  figure  for  the  record. 

(2)  I^Q  ^  T^o  and  I3Q  >  T3Q.  This  condition  probably 
indicates  noise  in  the  30-degree  channel.  The  data 
are  ignored. 

(3)  I^o  >  T10  and  ^30  —  t30’  This  condition  probably 
indicates  noise  in  the  10-degree  channel.  The  data 
are  ignored. 
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(4)  I10  >T^o  and  I3Q  >  T3Q.  This  condition  indicates 

detection  of  a  particle,  and  che  information  is 
stored  for  subsequent  treatment  along  with  other  such 
pairs  in  the  record. 

It  has  been  found  from  experience  that  the  occurrence  of 
conditions  (2)  and  (3)  is  relatively  infrequent  as  compared  with  the  occur¬ 
rence  of  condition  (1),  so  that  the  accuracy  of  establishing  a  background 
level  is  affected  little  by  the  fact  that  these  data  are  ignored.  However, 
it  should  be  pointed  out  that  this  is  a  consequence  of  selecting  the 
threshold  criterion  well  above  the  "noise*1  level.  Should  it  be  found 
desirable  in  the  future  to  reduce  the  threshold  criterion  (u.g.,  in  pur¬ 
suit  of  smaller  pulses),  then  this  question  should  be  re-examined. 

After  all  174A  data  pairs  of  the  record  of  interest  are 
examined  using  criteria  (1)  through  (4),  the  mean  background  values  for 
the  10  and  30-degree  channels  are  calculated  using  those  data  which 
satisfies  condition  (1)  above.  These  mean  values  are  then  subtracted 
from  those  raw  data  pair  values  which  satisfy  criterion  (4).  In  this 
manner  an  approximate  correction  for  background  level  is  applied  to  the 
pulse  height  data. 

The  data  pairs,  corrected  in  the  above  manner,  are  next 
sorted  with  respect  to  the  comparison  curves  (Figures  2  and  3)  in  order 
to  determine  the  particle  size.  Fir9t,  the  ratio  of  the  corrected  data 
pair  is  calculated,  and  compared  with  the  tabulated  values  which  were 
employed  to  generate  Figure  2.  If  the  calculated  ratio  is  greater  than 
the  value  corresponding  to  u  «  2  in  Figure  2,  it  is  assumed  that  the  parti¬ 
cle  is  small,  and  processing  by  the  absolute  intensity  method  is  initiated. 
However,  if  the  ratio  is  smaller  than  that  appropriate  for  a  *  2,  one 
initially  assumes  that,  the  particle  has  cs  >  2  and  the  ratio  method  is  appro¬ 
priate.  However,  before  it  is  finally  assumed  that  the  particle  is  large 
(o  >  2),  and  treated  appropriately,  a  check  is  made  to  insure  that  the 
amplitude  in  the  10-degree  channel  is  high  enough  to  correspond  to  a 
large  particle.  In  order'  to  satisfy  the  test,  the  amplitude  in  the  10- 
degvee  channel  must  exceed  the  threshold  criterion  by  an  order  of  magni¬ 
tude  which  has  been  Selected  somewhat  arbitrarily.  It  is  significantly 
less  than  the  expected  amplitude  in  the  10-degree  channel  for  particles 
with  a  »  2.0. 


If  the  amplitude  in  the  10-degree  channel  does  not  satisfy 
this  condition  (the  absolvite  amplitudes  increase  very  rapidly  for  particles 
of  0!  <2),  the  pulse  is  treated  using  the  absolute  amplitude  method,  thus 
implicitly  assuming  some  noise  in  the  30-degree  channel.  However,  if  the 
data  satisfy  both  the  ratio  and  absolute  criteria  for  larger  particles, 
the  size  is  determined  by  inverse  tabular  interpolation  using  Figure  2, 
and  the  value  is  stored  to  average  with  other  values  from  the  given 
record  and  other  records  being  averaged  for  the  desired  altitude  increment. 
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Assuming  the  the  corrected  data  for  a  given  pair  has  a 
ratio  greater  than  that  appropriate  for  a  =  2  particles,  or  if  the  ratio 
is  larger,  but  the  amplitude  in  the  10-degree  channel  does  not  exceed  10  T].q, 
the  size  of  the  particle  is  determined  using  inverse  tabular  interpolation 
in  Figure  3.  However,  data  treated  in  this  way  are  ignored  if  the  value 
of  a  determined  from  inverse  interpolation  from  Figure  3  exceeds  2.5.  No 
such  cases  have  been  encountered  so  far  in  the  data  analysis  performed 
to  date. 


A  schematic  flow  chart  may  help  in  a  detailed  understanding 
of  the  above  sequence  of  operations  given  in  Figures  4  and  5.  For  each 
record  processes,  a  count  of  the  method  ior  individual  pairs  is  maintained 
in  the  following  manner: 

(a)  Number  of  detected  particles  which  satisfy  the  ratio 
criterion  and  the  magnitude  criterion  for  the  10-degree 
channel  appropriate  for  large  (3  >  2)  particles 

(b)  number  of  detected  particles  which  satisfy  the  ratio 
criterion  for  small  particles,  and,  for  which  the 
assigned  a  is  less  than  2.0 

(c)  number  of  detected  particles  which  satisfy  the  ratio 
criterion  for  large  particles  but  for  which  the 
amplitude  in  the  10-degree  channel  indicates  that 
the  particle  is  small 

(d)  number  of  detected  particles,  which,  using  the 
absolute  intensity  analysis  method,  despite  ratios 
corresponding  to  large  3-values,  gives  values  of  Of 
exceeding  2.0  but  less  than  2.5.  For  these  detected 
particles,  the  bin  number  in  the  3  range  is  given 

so  one  can  see  what  3  value  is  assigned 

(g)  number  of  detected  particles,  which,  using  the 

absolute  intensity  method,  gives  values  for  a  ex¬ 
ceeding  2.  5. 

An  Investigation  of  the  data  records  processed  to  date 

ind icates: 


(1)  A  small  number  of  detected  particles  have  been 
treated  by  method  (a)  above. 

(2)  A  much  larger  number  of  detected  particles  have  been 
treated  according  to  p-escripcion  (b)  above. 

(3)  A  rather  large  number  of  detected  particles  have 
been  treated  by  prescription  <c)  above.  This  may 
indicate  calibration  or  threshold  errors. 

(4)  A  very  small  nuraber  of  detected  particles  have  been 
created  by  method  (d)  above. 

(5)  None  of  the  particles  have  been  treated  using  logical 
route  (o)  above,  to  date. 
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From  the  preceeding,  it  is  evident  that  the  data  have 
been  reduced  in  only  a  preliminary  fashion  to  date  and  that  additional 
numerical  experiments  must  be  performed  using  the  existing  flight  data 
in  order  to  optimize  the  procedure.  The  following  discussion  delineates 
several  specific  recommended  measures  in  this  regard. 

(i)  An  attempt  should  be  made  to  establish  a  pulse  criterion 
which  is  better  suited  for  the  identification  of  smaller 
pulses  than  that  presently  in  use.  Numerical  experi¬ 
ments  to  this  end  should  be  addressed  first  to  assessing 
the  sensitivity  of. the  derived  results  upon  this  param¬ 
eter. 

(ii)  The  data  reduction  method  should  be 

modified  in  order  to  better  allow  for  possible  dis¬ 
continuities  of  the  data  for  particles  near  a  =  2.0, 
where  reduction  procedures  change  from  absolute  to 
ratio. 

(iii) Use  of  the  30-degree  channel  in  place  of  the  10-degree 
channel,  for  the  absolute  size  determination  mode,  Bhould 
be  explored.  Selection  of  the  10-degree  channel  in  the 
present  scheme  was  made  on  the  basis  that  this  channel 
provided  greater  accuracy  in  the  determination  of  a, 

by  virtue  of  a  steeper  slope  for  the  output-versus- a 
curve.  However,  this  advantage  was  gained  at  the 
expense  of  greatly  reduced  sensitivity  for  smaller 
particles.  A  trade-off  in  the  reverse  direction 
might  well  be  desirable. 

(iv)  The  present  program  leads  to  incorrect  record  counts 
if  the  final  record  of  an  altitude  sub-average  for 
total  particle  concentration  shows  low  flow  rate 
(and  is  thus  ignored).  This  should  be  corrected. 

(v)  Modification  of  the  program  should  be  made  so  that 
the  total  number  of  particles  assigned  to  each  a- 
increment,  during  the  processing  of  those  records 
resulting  in  a  particle  size  distribution,  are  out- 
puted.  This  additional  output  datum  would  enable 
estimation  of  the  expected  errors  due  to  the 
counting  statistics. 

A  listing  of  the  currect  program,  together  with  a 
description  of  the  input  data,  is  given  in  Appendix  A. 
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C.  BALLOON  FLIGHT  RESULTS 


The  preliminary  results  of  the  28  October  balloon  flight  are 
presented  here  with  no  attempt  to  analyze  or  interpret  same.  Size  dis¬ 
tributions  for  3  to  10,  10  to  12,  12  to  14,  14  to  16,  16  to  18,  and  18 
to  20  kilometers  on  the  upleg  portion  are  presented  in  Figures  6  through 
11  respectively.  Figure  12  pertains  to  float  altitude  (i.e.,  >23  km), 
while  Figures  13  through  18  show  the  respective  downleg  size  distributions 
for  the  following  height  intervals:  22  to  21,  21.5  to  20.5,  20.5  to  19.5, 

19.5  to  18.5,  18.5  to  17.5,  and  17.5  to  16.5  kilometers.  Upleg  data  for 
other  portions  of  the  flight  have  not  been  included  here  owing  to  negative 
flow  rates,  as  discussed  previously.  In  addition,  owing  to  the  fact  that 
the  flight  lasted  longer  than  the  six-hour  recording  interval,  downleg 
data  below  16.5  km  were  not  recorded.  Note  that  all  figures  indicate  the 
presence  of  a  peak  at  Ct  =  1.1  -  1.2.  This  is  believed  to  be  a  consequence 
of  the  present  manner  in  which  threshold  criteria  are  applied  in  the 
computer  program,  as  discussed  in  the  previous  section  of  this  report. 

It  should  also  be  noted  that  a  tendency  is  shown  for  greater  spread  in 
the  data  for  that  portion  of  the  size  distribution  greater  than  about 
a  =  2.5.  Although  it  is  tempting  to  ascribe  this  to  statistics  on  the 
ground  that  smaller  total  counts  have  been  recorded  in  this  size  regime, 
the  transition  is  too  abrupt  to  be  entirely  plausible. 

Figure  19  is  a  plot  of  total  particle  (i.e.,  between  0.2  and 
1  micron  diam)  concentration  profile  observed  during  this  flight.  This 
shows  the  well-known  peak  at  about  20  km  altitude,  as  well  as  the  trough 
between  12  and  18  kilometers,  together  with  the  rapid  fall-off  with 
increasing  altitude  below. 

The  lack  of  overlap  between  upleg  and  downleg  data  is  an 
unfortunate  consequence  of  the  low  flow  rate  conditions,  which  necessi¬ 
tated  the  rejection  of  certain  entire  portions  of  the  total  concentration 
computation.  The  effect  on  size  distributions  determination  on  the  other 
hand,  is  not  as  severe.  These  data  (namely,  16  to  18  and  18  to  20  kilometer 
data)  have  been  retained,  although  the  reader  is  cautioned  as  to  the  possible 
unreliability  of  the  ordinates  in  these  two  instances. 
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Figure  19.  Altitude  profile  of  particle  concentration  within  diameter 
size  range  0.2  -  1.0  micron. 
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APPENDIX  A 
PROGRAM  LISTING 


REAL  IH0USE*IDATA»ITRESH»JT3E5H 
INTEGER*2  JDATA 

DIMENSION  HOUSE (32) ♦ I  DATA (3488) *AVGlO<91) »AVG30(91) 
DIMENSION  RATI0(9I) »ST0RE(9i) *DATA(3488) »BIN(9D 


DlMtNSION 

dimension 

dimension 

DlMtNSION 

COMMON 


JDATA(1760) 

<DAT  A ( 1 760 ) 

DAR (3488) 

NSUB (9) 

z  *flowrt  *alt*alaser*photoi *PH0T02  *HV  *BATVOL*BATAMP « 


ltLECT£*HLA$TE*A:_A5C*VOLTL*VOLTH 


sroKt.2=o.o 

DO  11  1=1*91 
STORE ( I ) =0*0 

n  continue 

IRtCN0=0 
RE w I NQ  15 

C  READ  AVG10  AND  AV630 

READ ( 5 ♦ 100 )  RC0N5T * ACONST 
100  FORMAT  (2F10.0) 

READ (5 ♦ 10)  AVG1C 
10  FORMAT  (5E14.5) 

READ (5*10)  AVG30 
DO  a  1-1*91 

RAT  10(1)  =  (AVG30(I)/AVG10(I) >*RCDNST 
AVG10  ( I )  sAv/GIO(  I )  *  A  CON  ST 
71  CONflNuE 

RkITE(6*lO)  AVG10 
ARITfc(b*10)  AVG30 
aRI  ft  (6*910)  RATIO 
910  FORMAT  U0E11.3) 

C  READ  the  NUMBER  of  DATA  BLOCKS  TO  BE  PROCESSED  AND  THE  TRESHOLDS 
READ (5*4621)  IREC* ITRESR * JTRE5R 
4621  FORMAT (515) 

60  FORMAT  (I5.2F5.0) 

C  ROAD  THE  HOUSEKEEPING  PORTION  OF  THE  DATA  BLOCK  AND  THE  DUMMY  LEADER 

2700  lNDi=0 

2701  READ ( 5  *2702 )  I T1 *  I  ACT 
READ ( 5 *2 702 )  NSJBI 

RE AD ( 5  * 66Q0 )  ( NSUB ( I ) *I  =  1»NSUBI ) 

6600  FORMAT (1014) 

C  NSUB  IS  HE  NUMBER  OF  RECORDS  TO  BE  PROCESSED  FOR  FINER  ALTITUDE 
C  TREATMtNT  FUR  TOTAt  PARTICLE  NUMBER 
1565=1 
ISU98=0 
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iREClN=IRECNO+l 
IR£CFN=IRECIN+N5UB( 1565) 

198  =  NSUBU565) 

1570=0 

3SJM=0.Q 

IPRT=0 

2702  FORMAT <2l4) 

2703-  CONTINUE 

READ! 15*700»END=198*ERR=112)  JDAT  A 
IRECNQ-IRECNO+1 
700  FORMAT (6 (255A2) »23CA2) 

IN01SIND1+1 
1570=1570+1 
2003  FORMAT ( 1H  *20165 

IF(IACT)2704»27o5«2705 
2704  IF (INDI-IT1) 2703 12700*2700 
27J5  CONTINUE 


*R  HE  (6*2912)  (JOATAU)  *1  =  1*48) 

2912  FORMAT (8110) 

00  3000  1=1*1760 

IF ( JDAT A ( 1 ) ) 3001 i 3002 i 3002  ..  .  _  _ 

3002  4DAT4(I)=J3ATAU) 

00  10  3000 

3001  <DATA(I)=J0ATA(I)+65536 
3000  CONIINUE 
00  10  141 
198  WRITt (6*2004) 

2004  F ORMAT ( 24H  END  OF  FILE  *  IRECNO  =  ) 
wRI  TE (6*2003)  IRECNO 
1RECN0=IRECN0+1 
1570=15  ’0+1 
IN01  =  IND*"  I 

IF  UACT)  270/ >2708*2708 

2707  IF ( INUl-IT 1 )  27U3 1 2 700  * 2700 

2708  IF(IMD1“IT1)  2703«2706*2706 
112  IR&CNO= IRECNQ+1 

WRITE (6* 1004) 

1004  FORMAT (56H  A  TAPE  ERROR  IN  RECORD  HAS  BEEN  DETECTED*  DATA  IGNORED 
1) 

IN01MND1  +  1 
1570=1570+1 

I F ( I  ACT  >  2 709*2710*2710 

2709  IF(lNDl-ITl)  2703*2700*2700 
2 7 1Q_  IFUNDl-ITl)  2703*2706*2706 
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141  J=1 

DO  151  1=1*16 
IH=<DATA  (I  )  /  2  56 
IHOOSE ( J)  =  h 
IH=KDaTA (I ) -IH*256 
iHOjbE ( J+l ) s|H 
J  =  J*2 

151  CONTINUE 
J  =  i 

DO  161  1=17*1760 
ID=<DATA ( I ) /256 
1DAT  A ( J) = I D 
ID=KDATA(I ) -10*256 
IDATA(J*1)=ID 
J  =  J*  2 

161.  CONTINUE 
/  1005  FORMAT (6 (£18*4) } 

CAuL  EVALJA(IRECNO*lHOUSe> 

I F (I  HOUSE (17) -46 • )  7500*7501*7501 
T500  1*6=198-1 

IF ( INDi-ITl )  2703*7502*7502 
7502  GO  TO  47 
7501  CONTINUE 

C  SCALE  THRESHOLDS  TO  ALLOW  FOR  LASED  POWER 
FACTORS  I  HOUSE (19) /FLOAT (102) 

Y10=FL0AT(ITR£SR) 

Y3U=FL0AT(JTRES9) 

X10T  =  C0NVRT  < Y 10) ^FACTOR 
X30T =CCNVRT (Y30) *FACTOR 
YlO=(7.984375+ALOG10(XiOT) )/0.015625 
Y30=( 7.984375* ALOG10(X30T) )/0. 015625 
I TRfcSH=Y 10 
JTRtSH=Y30 

WRI TE (6*4008)  ITPE5H* JTRESH 

4008  FORMAT (20H  NEw  10  THRESHOLD  =  *F10.4*20H  NEW  30  THRESHOLD  =  *F1Q.< 


1) 

44  ALTITD=ALT 
114  <0 JNT  s 1 
I  I o\  =  0 
'  KTR  =  1 
<0JNTR=0  • 

DAT  AR 1 =0*0 
DATAR3=0*0 
aR 1 TE ( 6  *920) 

920  FORMAT  (34H  DATA  HAS  BEEN  READ  FROM  THE  DISK) 
DO  51  1=1*91 
BIN(I) =0.0 
51  CONTINUE 
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.v.\' J.V>  ii IFA  a.V  &*&*»& 


C  CHECK  FOR  PARTICLES 
JO  2 1  1=1*3488*2 
IF  (iDATA(I) -I TRESH)  4000*2  3*23 
23  IF  (IDATA( 1*1) -JTRESH) 4002 *24*24 
4Q00  IF(IDATA(I*1)-JTRESh)  22*22*4001 

4001  I  l3fN=I  IGN+I 
-00  TO  21 

4002  IF(iDATA(I)-lTRE5H)  22*22*4003 

4003  1 1 GN  = I IGN+1 
—GO  TO  21 

2  4- I  F  ( 1  DAT  A  { I )  -256* )  25  ♦2  0»2,O 
25-1 F (lDATA(I+l)-256*)  26*2«*2€ 

26  DATA(K0UNT)=C0NVRT(IDATA(I) ' 

DATA (K0UUT*1)=C0NVRT<IDATA( 1*1) ) 

K0JNT=K0UNT*2 
-GO  TO  2$ 

22-IFUDATA(I).EQ.0.0)  GO  TO  20 
“I F ( 1  DATA ( I ♦ 1 ) .EO.O.C)  GO  TO  20 
OAR  (<TR)  =CDMVRT  (IDATAU  )  ) 

DAR(<TR*l)  =  CONVRT ( IDATA ( I *1 ) ) 

QATARl=DAR (<TR) +DATAR; 

DATAR3  =  DAR  (<TW  +  U  >DATAR3 
<TR=<TR*2 
<0JNTR=K0UNTR*1 
21  CONTINUE 

1 F ( <OuN  TR )  2000*2000*2001 

2000  DAT AR i =0*0 
DA  l AR j=0.0 
RAT  1 02  =  0*0 
GO  TO  2002 

2001  DATARi^QATARl/OUNTR 
DATAR3=DATAR3/K0UNTR 
RAT  102  =  0*0 

2002  CONTINUE 

WRITE  (6*4004)  DATAR1 *0ATAR3 

4004  FORMAT ( 1 7H  MEAN  10  3KGND  =  *El5.8»\7rt  MEAN  3u  BKGND  =  «E15.8) 

KOjMT=KOjNT-2 

<TR=<TR-2 

<<<T=K0UNT/2 

/(RITE  (6.4Q05)  <<KT  .KOUNTR  « IIGN 

4005  F  QR  -1 A  T  (28h  ^ARTICLFS  IN  THIS  RECORD  =  *I4«15H  BKGND  COUNT  =  *14. 
1 1  7R  IGNORED  PAIRS  =  *14) 

C  CAlCUuAT£  STANDARD  DEVIATION  DF  BACKGROUND 

SUM  1=0*0 

6UMt=0.0 

•  DO  h215  I  =  1  *<TR  *2 
SUM1 =SUM1 ♦ ( DAK ( I ) -QATAR  1) **2 
SUM2  =  S'JM2 ♦  (  DAR  ( I  ♦  1 )  -QATAR 3)  *«2 
4215  CONI iNJt 


55 


5rA\Pl  =  i>URT  (5UM1/K0JNTR) 

STAND3=SQRT(SUM2/KOJNTR) 

XiO=lTRESH 

X3Q=JTft££ri.  _ _ _  _  _ _ _ _  . 

xi  hc=cqnvrt  (xio) 

XJTHC=CUNVRT (X33) 

WR I T  E ( 6  ♦ 4006 )  XiTHC.XJTHC 

4006  FORMAT (26H  CI.NVERTED  10  THRESHOLD  =  «E15#8.26H  CONVERTED  30  THRESH 
lULD  =  « c  1 5  •  8 ) 

WRITE (6*4007)  STANDI .STAND3 

*♦007  FORMAT  ( 20H  10  0<GN  STND  DEV  =  »£15.8«20H  30  BKGN  STND  DEV  =  .E15.8 
1) 

C  CALCJLATE  AND  APPLY  RAYLEIGH  SCATTERING  CORRECTION 
WR I TE { 6  *940 ) 

940  FORMAT  ( 32H  FINISHED  CHECKING  FOR  PARTICLES) 

ICR  =  0 
I C  A  =  0 
ICRa=q 
ICAR=0 
ICARUO 

DO  31  I  =  1  ♦  <  3JNT  »2 
I FLoG= ) 

OATA ( I ) =DA  T  A ( I ) -DATAR1 
DAIA(  I*1)=DATaU+D-DATAR3 
C  CAt.CJLATE  THE  RATIO  AND  SORT 
ARAT IQsQATA ( I ♦ 1 ) /DATA! I ) 

DO  61  J= 1 ♦ 9 1 

IF (ARaTIO-RATIO(J) ) 61 » 36 » 36 
36  IF (J-30> 32.32*4018 
32  CALL  HECALOATA(I)  .AVG10) 

DO  171  IJ=1.91 

IF (DATA (I ' -AVGlOt I J) )  1 72.1 72.171 
1 t 2  J=I J-l 

IF (J) 31.31.315 
315  IFU-30)  *.009.4009,416 
<.009  IF  (  1F[  oL*  1  )  40  1  1  .AO 1 0.40 II 
**013  ICAsICAM 
GO  TO  37 

4011  ICRAs  J  CRA*1 
GO  TO  37 

416  IFU-40)  4012. 4012.4013 

4012  WRITE (6.4014) j 

4014  FORMAT (3  3m  OVERLAP.  PARTICLE  COJNTED.  J  =  .14) 

IF(IFlGG-I)  4015.4016.4015 

4016  I CAH=  J  CAR ♦ 1 
GO  TO  37 

4015  ICRA=ICRA*1 
GO  TO  37 

4013  WHI Tfc (6.4017) J 

4017  FORMAT ( 26H  toRG  OVERLAP.  IGNORE*  J  =  *14) 

ICAR  1 s I  CAR  1 *1 

GO  TO  31 


n  n  r>  rv 


1 n  CONTINUE 

4018  CONTINUE 

C  bE  IF  AMPLITUDE  IS  LARGE  ENDJGH 
5ECCT=XITHC*10.0 
IF (OAT A! I ) -SECCT) 4019*4019,4020 

4019  I FlGG=0 
GO  TO  32 

4020  ICR=ICR+1 

go  ro  37 

37  BIN(J)  =  BIN(  J)  n*0 


GO  TO  31 
61  CONTINUE 

RAT 102=RATI 32*1.0 
31  CONTINUE 

WRITE (6,950) 

960  FORMAT  (20H  ALL  DATA  lb  SORTED) 
wKi Tt (6*4021) 

4021  FORMAT <  3 IH  FOR  THIS  RECORD  THE  COUNTS  ARE) 

WRITE (6*4022) ICR 

4022  FORMAT  (3*»H  ^ARTICLES  ^ROM  RATIO  TEST  ONLY  =  *14) 
wRl  IE (6*4023)  ICA 

4023  FuRRAT ( 33H  PARTICLES  FROM  aBS.  CAL.  ONLY  s  ,14) 
rtRl TE (6,4024)  ICAR 

4024  FORMAT (38H  PARTICLES  MOVED  FROM  ABS.  TO  RATIO  =  ,14) 

WRITE (6,4025) ICRA 

4025  format ( 3hh  particles  moved  from  ratio  to  abs.  =  *i4) 

wRl IE (6,4026)  ICAR  I 

4026  FORMATOOH  TOO  LARGE  OVERLAPS ,  I  SNORED  =  ,14) 

CORRECT  THE  NUMBER  UF  COUNTS  IN  EACH  B.’N  DUE  TO  INSTRUMENT 
DEAD  TIME  AND  FLOW  RATE 

CORRECTION  DUE  TO  INSTRUMENT  DEAD  TIME  IS 
lRRr=lCRMC4*ICAR*ICRA 
COJNT  =*U'JNT/2 
COUNT  - 1 744. -COUNT 

COR  I  NS=  (U6./ll.)»l  744.*AL0G(  1744. /COJNT)  /  ().  744  .-COUNT  )  ) 

US JMsQS JM.F^OAT ( IPRT) / ( F lOWRT* 1 7. 44) *COR INS 
IF ( 1 570-NSU3 ( I  565) ) 6500,6  501,6501 
6501  I F (  1 9  8 )  6790,6791,6  790 
6  791  USU'i  =  0»0 

GO  TO  6792 

t  790  JSJM=OSJM/F„OAT< 198) 

6792  CONTINUE 

wHl IE (6,6502)  IRECIN,IRECFN*198 

WRITE (6,3581) OS JM 

USJM=0.0 
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IPRf=0.0 
1565=1  565  +  1 
lb  70*0 

ISJ98=ISU98.I98 

IRtClN=IRECN0.1 

1RECFN=IRECIN*NSUB(I565) 

65q2  FORMAT (22H  INITIAL  RECORD  NO.  =  »I4»2qH  FINAL  RECORD  NO.  =  ,14,20 
l  RECORDS  INCLUDED  =  »I4) 
bbOO  CONTINUE 

C  CORRECTION  DJE  TO  FLOW  RATE 
COREL 0=1*/ (FlOWRT*I7*44) 

1)0  41  1  =  1*91 

STORE ( I )  =8IN(I)*C0«1NS*C0RFL0+ST0RE(I) 

41  continue 

ST3RE2=STORE2.RATI02*CORINS*CORFLO 
IF  (  IND1-IT1 ) 2703,2706,2706 
270b  CONTINUE 


121  wRITE(6«930)  IRF.CNO 
930  FORMAT  (HO) 

47  CONTINUE 
43  WR 1 TE { 6  »40 ) 

40  FORMAT  (i*l**/V/»40X.*THE  DISTRIBUTION  OF  PARTICLE  SIZE  IS. *///♦ 
ISO*  A„PhA  NJMBER  OF*) /t5(15X»*PARTICLES*)  ) 

ALPHA=  ,45 
00  131  IJ-1.91 
ALPhAs ALPHA. ,05 


IF ( 15j9d) 6794,6793.6794 

6793  STORE ( I J) =0.0 
SO  TU  6795 

6794  STORE  (I  J)  sSTOREUJ) /FL0ATUSU98) 

6  79 5  CON  1 iNUt 

*RITE(6»50)  AlPHA.STORE ( I J) 

SO  FORMAT  !5 (F12.3.E12.3)  ,/) 

131  CONTINUE 

IFUSU98)  6796.6797,6796 
6797  i>T0RE2  =  0.0 
CO  TO  6  796 

6796  3TUR£2=STORE2/FLOAT<ISU98) 

6796  CONI INUt 

WRITt(6»90)  ST0RE2 

90  FORMAT  (/,*  ThE  NUMBER  OF  PARTICLES  DETECTED  OF  DIAMETER  GREATER  t 
IRAN  1  MICRON  IS.*E14.3) 

3581  FORMAT  ( 29.H  MEAN  NO.  OF  PART1CLES/CC  IS  ,F8.4) 
lTII JysALTI TO 
00  3571  1=1*91 
STORE (I >=IOOO*5TORE<!) 

35  71  CONTINUE 

CALL  HISTURECN0.ST0RE.90) 
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o  o  f  i  o  o  o  o  n  r>  o  r^n  n  o  o 


DO  3b 72  1  =  1*91 
IF ( STORE ( I ) )  3572*3572*3573 
3573  CONI 1NUE 

STORE ( I) =AL0G( 10 . 0*STORE ( I) ) 

3572  CONTINUE 

CALL  HIST (IRECN0»ST0RE*9Q) 

00  111  1=1*91 
STOWE ( I ) =0.0 
111  CONTINUE 
STOWE2=0.0 
GO  TO  2700 
45  CONI INUE 
STOW 
END 

SUBROUTINE  EVALJA(IRECNO«lHDUSE) 

THIS  PROGRAM  C0NVERT5  THE  H0USE<EEP1NG  DATA  INTO  THEIR  PROPER  UNITS 

WEAL  ihouse 
DIMENSION  I  HOUSE ( 32) 

common  2*flowrt*alt*alaser*photoi«photo2«hv»hatvol*batamp* 

1ELECTE  *ALAST£*ALASC»VOLTt*V0LTH 
THE  VARIABLES  REPRESENT  THE  FOLLOWING 


8ATV0L  =  BATTERY  VOLTAGE 
FLOaRT  =  Fv-Ow  RATE  BATAMP  =  BATTERY  CURRENT 

Alt  S  ALTITJDE  ELECTE  =  ELECTRONICS  temp 

ALAScR  =  LASER  POWER  ALASTE  =  LASER  TEMP 

PHOTOl  =  TEMP  OF  PHOTOMULTIPLIER  NO  1  ALASC  S  laser  current 

WHOT02  =  TEMP  OF  PHOTOMULTIPLIER  NO  2  VOLTL  =  OV  REFERENCE 

HV  =  HIGH  voltage  volth  =  sv  REFERENCE 

TEMP  IS  IN  DEGREES  AflSQLUlE 
wWITF(6»b00>  IHOUSE 
600  FORMAT  Ht20*3) 

IF ( IHOUSt (20) ) 5*5*6 

6  L- U 2.04*2. 861* (5. * I  HOUSE (20) 7256.)  ) 

5  IF (lHOJSt (17) ) 15,15*16 
16  FV  =  S.0«lH0JSEU7)/256* 


F V2  =  F  y«F  V 

FLO, iRT  =  -64, 39*137, 23 *FV-5»666*FV2 
Fl0awt=FL0WRT/6O. 

15  IF ( HOUSE ( IB) .E0.O.O)  GO  TO  25 
ALT=  (o*  10*5.*  HOUSE  (  IB  >  7256.) 

25  IF  UhnuSE  (  19)  .EO.O.u)  GO  TO  35 
MLASERslHOJSE (19)»30,/256. 

35  IF  U  house (21) ) 45 ,45  *46 
«*b  PH01  Ol  =  TH  (  HOUSE  (21)  ) 

*♦5  IF  1 1  HjuSt  (22)  >55*55*56 
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O  C\ 


bfe  PHD  I Oc  =  TR ( I  MOUSE  (22 )  ) 
t)5  HV=25oO.*IH0JSE(23)/256. 

BATv0L=-50.*IH0U5E(24)  /256. 

BAT«Mp=io,#IHOU5E(25) /256. 

IF  (  HOUSE  (26)  ) 65.65.66 
66  tlECTE  =  TK  (  HOUSE  (26)  ) 

65  IF  (  HOUSE  (27)  )  75 . 75 .76 
fb  AL«STE  =  TK  (  HOUSE  (27)  ) 

75  AlASC  =  330.*I HOUSE (28) /2 5 6 • 

VOL  TL  =  5 • *  I  HOUSE (29) /256 • 

V3LH=5»Q40#I HOUSE ( 30 ) /256. 

14  *RITL(6.30) 

30  FO^AT  (  HI  «///.20X»*'ThE  HOUSEKEEPING  DATA  HAS  THE  FOLLOWING  VALUE 


WR  1  ?  E  ( 6  *  40)  Z.FlOWRT«ALT.ALASER 

40,J,?2viTy::  '"-miTER  '88.3,  8X.XFL0.,  RATERS. ALTITJ 

lDE*‘Fb.3*i9X*^LASER  POWER^FB.3)  ^LIITU 

«R I t t ( 6 *50 >  photo i ♦ photo 2 *hv*batvol 

5°  RORv,AT  (/,*  PmOTOTUBE  1  TEMP.  *F8  .  3  *  10X  ,  ^PHOTOTUBE 
1  Hio)  VCLTAGL»lfF8.3*15XtJ'BATTF.RY  VOL  TAGE  8 • 3 ) 

1  T  E  (6*60)  BATAMP.ELECTE.ALASTE.ALASC 
so  FORMAT  </.*  BATTERY  CURRENTS.  3. 12X.HELECTRUNICS 
1  TfcMPERATURE*F8.3»10X. GLASER  CURRENT AF&.3) 

1 T  E ( 6 . 7o  VOLTL»VOLTH 

'DjFO-HaT  </,*  ov  REFERENCE  VOlTAGE*F  8.3.7X.X5V  REFERENCE  V0lTAGE*F8. 
13  Re  )  v'R'M 


2  TEMP*  *F8 •  3 » / 


TEMPi*Fa.3./« 


tND 

Function  cdnvrt<  a> 

C0NVRT= (10.**(-(255.-A)/64,> )*i.0£-Q4 
return 


♦12. 7R  POINTS./) 


END 

SUBROUTINE  RI5T (NU.FRE0.IN) 

DIMENSION  UDJT  UOO)  .FREQUOO) 

DATA  K / IH*/ .NOT R/ H  / 

1  FORMAT  (6N  EACH  »A1.8H  EQUALS 

2  FORMAT  (  J&.4X.  VODAI) 

3  FORMAT  (  NhQ  I  NTt’S/A^  *19(12.3X1 
h  FORMAT (1H1.A7X.ZIM  FTNAL  RECORD 

5  FOR  RAT  1 1CH0FREQUENCY) 

6  F OR hat  (6h  CLASS) 

7  <■  <113-  . - 1- - - - 

10  FORMAT  (  H* )  ** 


.12) 

NO,  IS 


H) 


P«1NT  TITLE  AND  FREQUENCY  VECTOR 

..RITE  (6,4)  NJ 
DO  12  I  =  I,  I N 
i?  JOJT  (  !  ) =FR£0( I ) 


SO 


nrm  n  on  non  o  n  o  nnn  non 


rtRl TE (6.5} 
wRl TE (6,7) 


FIND  LARGEST  FREQUENCY 


FMAX=Q.0 
DC  20  1  =  1,  IN 
IF (FREQ( I ) -FMAX) 20,20,15 
15  rMAX=FREQ(I ) 

20  CONTINUE 


SCALE  IF  NECESSARY 


JSCAL=I 

IF (FMAX-50. 3) 40,40,30 
30  jSCALs (FMAX+49.0) /50.0 
ttR I T  E ( 6  » I )  < .JSCAL 

CLEAR  OUTPJT  AREA  TO  0LAN<S 

40  DO  50  1=1. IN 
50  JOJI (I)sNOTR 

LOCATE  FRFOJE'NCIES  in  EACH  interval 

MAX =F MAX /FLOAT ( JSCAL) 

DO  bO  1=1. MAX 
X  =  M  A  X  — ( 1-1 ) 

DO  70  0=1 .  IN 

1F(MEq{j)/fIQAT{jSCAL)-X)70.60»6Q 
60  JOJI ( J) -K 
TO  CONTINUE 

IX=X»FLOAT ( JSC  A; ) 


PRINT  LINE  OF  FREQUENCIES 
80  wR l T  c (6.2  )  IX.(JOJTIJ) .J=1.IN) 

generate  constants 

jo j  t { n  =  1 
Js  i 

DO  90  I  =  5,  IN,  5 
J  =  J*1 

90  jOJT ( J) -  I 
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PRINT  INTERVAL  NUMBERS 


wP  !  TE  {6  t  7 ) 

WRITE (6  1 3 }  ( jOUT { J  5 «Jsl»19) 
wR I T  E ( 6  1 6 ) 

RETJRN 

ENP 

SUBROUTINE  RECAL {DATAl , AVGIO) 

DIMENSION  AVGIO (913 

COMMON  Z  ♦FLOa'RT  %  aLT  »  Al  ASER  »PHDT01  ♦  RMDT02  »Hv  *8ATVOL .BATAMP  , 

1  ELEC TE*AL AS TE  tALASC* VOLT L*VOLTH 

DATAl=DATAU5.F  +  6»2.33/{ALA5ER/6,> 

RETURN 

END 

FUNCTION  TR(Ab) 

T=ALOG( (256. -AS) /AB) 

TR=298./ ( I.+.0BO6»T)-273. 
return 


END 
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